Replica theory of the rigidity of structural glasses 
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We present a first principle scheme to compute the rigidity, i. e. the shear-modulus of structural 
glasses at finite temperatures using the cloned liquid theory, which combines the replica theory and 
the liquid theory. With the aid of the replica method which enables disentanglement of thermal fluc- 
tuations in liquids into intra-state and inter-state fluctuations, we extract the rigidity of metastable 
amorphous solid states in the supercooled liquid and glass phases. The result can be understood 
intuitively without replicas. As a test case, we apply the scheme to the supercooled and glassy state 
of a binary mixture of soft-spheres. The result compares well with the shear-modulus obtained by 
f^ ■ a previous molecular dynamic simulation. The rigidity of metastable states is significantly reduced 

with respect to the instantaneous rigidity, namely the Born term, due to non-affine responses caused 
by displacements of particles inside cages at all temperatures down to T = 0. It becomes nearly 
independent of temperature below the Kauzmann temperature Tk. At higher temperatures in the 
Pj ' supercooled liquid state, the non-affine correction to the rigidity becomes stronger suggesting melt- 

^ I ing of the metastable solid state. Inter-state part of the static response implies jerky, intermittent 

stress-strain curves with static analogue of yielding at mesoscopic scales. 
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PACS numbers: 61.43.Fs 



I. INTRODUCTION 

Rigidity is a distinct character of solids including crystals and glasses [l|. It is quantified unambiguously by the 
shear-modulus [2| which represents strength of materials against shear deformation, which changes the shape of the 
containers but not their volume. The shear-modulus is finite in solids but zero in gasscs and liquids so that, much 
like the order parameters, it distinguishes solids from other states of matters, in sharp contrast to the bulk modulus 
which is finite even in gasses. We denote the shear-modulus simply as rigidity in the present paper. It is the rigidity 
of a Goldstone mode which reflects the breaking of the translational symmetry much as the rigidity of the spin- wave 
in ferromagnets which reflects the breaking of the rotational symmetry. The significance of the shear-modulus can 
be hardly appreciated at the level of the macroscopic elastic theory which treats all clastic constants on the same 
footing. We must go back to the microscopic scales to study the mechanism of the rigidity of solids. Unfortunately the 
rigidity of solids is usually taken for granted and regarded merely as a fitting or adjustable parameter. In the present 
paper we dig into this basic problem and attempt to develop a theoretical approach to compute the shear-modulus 
of structural glasses at finite temperatures starting from microscopic Hamiltonians. Microscopic computation of the 
rigidity of structural glasses would be useful also for practical applications of amorphous materials. 

A class of amorphous solid or glasses can be obtained through supercooling liquids. First one cools a liquid down to 
the supercooled liquid state, the metastable liquid phase below the melting transition temperature T^. A supercooled 
liquid behaves as a visco-elastic material. On one hand the system behaves as a solid with a finite effective rigidity 
for a long time before it flows. Within this time scale, which is called as the /3-regime, the configuration of the 
particles does not evolve much but the particles fluctuate mostly within narrow spaces around themselves called as 
cage created by the surrounding other particles. On the other hand the system behaves as a fluid with high viscosity 
at longer time scales called as the a-regime. In this regime the configuration of the cages themselves arc reorganized. 
The two qualitatively different relaxation mechanisms manifest themselves as two step relaxations, i. e. a and^- 
relaxations, in various observables such as the intermediate scattering function, dielectric susceptibility and e.t.c. [3|. 
Upon lowering the temperature, the viscosity or the structural relaxation time t^ increases enormously. Eventually 
Ta exceeds the laboratory time scale at some glass transition temperature Tg and the system falls out-of equilibrium. 
As the result we are left with a piece of glass with a finite rigidity, which behaves as a solid in all practical means. 
This is the glass transition observed in practice. Thus the rigidity of glasses is completed by enormous increase of the 
viscosity in sharp contrast to crystals where the rigidity is established abruptly by 1st order phase transitions with 
no appreciable precursors. The rich visco-elastic or rheological properties, which are absent in crystals, give distinct 
flavor to supercooled liquids and glasses and make them useful in practical applications. Apart from the interests to 
the outstanding question that whether an ideal, thermodynamic glass transition exists or not, development of flrst 
principle, microscopic theories to describe such an out-of equilibrium amorphous state of matters is strongly desired. 

Among various theoretical proposals Q, the Random First Order Transition (RFOT) theory p-Q provides a useful 
working ground to study the supercooled liquids and amorphous solids in a uniflcd manner. The basic scenario 
of the RFOT theory goes as follows at the mean fleld level. By lowing the temperature down to the so called 
dynamical transition temperature Td, exponentially larger number of metastable amorphous solid states emerge so 
that the structural entropy or the complexity becomes flnite. The ideal, equilibrium glass transition takes place at 
the putative Kauzmann temperature Tk due to the entropy crisis mechanism |1C| . Within the RFOT theory, T^ is 
equivalent to the critical temperature Tc predicted by the mode coupling theory (MCT) where Tq, diverges [ll|, [l^] ■ 
As is well known, the dynamical transition itself happens only at the mean-fleld level and T^ would remain in real 
systems at most as a crossover temperature. The metastable solid states can have only flnite life times in flnite 
dimensional systems. True divergence of the viscosity, if any, may occur approaching Tk as envisioned by Adam, 
Gibbs and DiMarzio[lJ,[lJ]. An important concern of the RFOT theory is to seek for a theoretical foundation of the 
scenario 0, i, [il-[i3 ■ 

At the mean-fleld level, the RFOT theory is backed up by the cloned liquid method [18l - l2l| and the modc-couping 
theory (MCT) [ll|, 113 • Consistency between the two approaches is currently investigated intensively [23] . The cloned 
liquid method provides a first principle scheme to compute the equilibrium and guasi-equilibrium properties of the 
supercooled liquids and amorphous solids starting from microscopic Hamiltonians with the help of the liquid theory 
[23| combined with the replica method [2J,l2a]- In the nutshell the cloned liquid theory, in its simplest formulations, 
views the metastable amorphous solids through an effective Einstein model in which each particle is subjected to a 
virtual Hookian spring with one end fixed at the mean position of each particle. The configuration of the latter is 
just that of a liquid. The virtual Hookian spring is meant to mimic the cage effect: the spring constant is inversely 
proportional to the size of the cage A, which is the order parameter of the theory. Physically, existence of amorphous 
states with finite cage size A means breaking of the translational symmetry. 

In the present paper we study rigidity of structural glasses against shear at finite temperatures using the fluctuation 
formula of the rigidity |26| . We evaluate it using the cloned liquid theory. This would amount to a formulation of an 
effective Dcbye model starting from the effective Einstein model, which is a necessary step to go beyond the mean-field 



theory. Long wave length elastic deformations are the essential low energy excitations in solids and play important 
roles much as spin- waves in ferromagnets. For example localized plastic events which are elementary steps of the flow 
in glasses |27H29l | are known to accompany smooth elastic deformations around them, which naturally implies elastic 
free-energy barriers and interactions with each other through the long-ranged elastic couplings |30l433| . 

We show that the cloned liquid method enables decomposition of a generic response function into two parts: 
the intra-state responses corresponding to the /3-relaxation and inter-state responses, which correspond to the a- 
relaxation. Based on this scheme we compute the intra-state rigidity /i, which presumably represents the effective 
rigidity of supercooled liquids and the rigidity of glasses. On the other hand, the inter-state fluctuation reduces the 
rigidity down to 0. In general, it is reasonable to expect that the /3-relaxation only weakly changes between ideal 
equilibrium and out-of equilibrium situations while the a-relaxation strongly changes j34| . Thus the quasi-static 
approximation for the intra-state responses would be valid for both the equilibrated supercooled liquids and out-of 
equilibrium glasses. 

We find that the intra-state rigidity jj,, i. e. the rigidity of the metastable glassy states is significantly reduced with 
respect to the Born term j35| which describes the instantaneous, affine response by a non-affine correction term which 
represents stress relaxations caused by relaxation of particles inside cages, i. e. the /3-relaxation, at all temperatures 
down to the zero temperature limit. Our result implies a characteristic temperature dependence of the intra-state 
rigidity. On one hand, it is nearly constant below the Kauzmann temperature Tk- On the other hand, it strongly 
depends on the temperature above Tk. To test our scheme, we applied the method to a binary mixture of soft spheres. 
We found the theoretical prediction compares well with the previous result by Barrat et. al. (1988) [Sa] obtained by 
a molecular dynamic simulation. 

An interesting general problem is how glasses melt. It has been proposed a long time ago by Born |37| that melting 
of solids may be signaled by vanishing of the rigidity. Although this rigidity crisis scenario obviously does not apply to 
the equilibrium liquid-solid transitions which are 1st order phase transitions, whether it is relevant for the melting of 
superheated metastable crystals approaching the spinodal temperature Tg from below is an intriguing question |38l - t4ll | . 
Interestingly enough, there is an intimate analogy [421 between the melting of metastable amorphous solids at the 
dynamical transition temperature T^ (or the MCT critical temperature Tc) and the melting of superheated metastable 
crystals at Tg. At the mean-field level T^ is regarded precisely as the temperature above which the metastable amor- 
phous solids states become absent in equilibrium. Moreover it has been pointed out that the underlying mechanism of 
the dynamical critical phenomena found by MCT is the qualitative change of the free-energy landscape at around Tj 
[8, 43]. Then an interesting question is whether the rigidity crisis scenario also applies to the melting of the amorphous 
solids. 

For clarity let us note that there are important classes of amorphous solids other than the glasses obtained by 
supercooling simple liquids. In a class of systems including gels, polymer glasses and rubbers, formation of disordered 
networks of molecules and colloids are important. Another important class is the granular matters which are athermal, 
i. e. the temperature play no essential roles. How the rigidity emerge in these amorphous systems have also attracted 
a lot of interests and some microscopic theoretical approaches have been developed [4J-[4g| . 

The organization of the paper is as follows. In the next section, we discuss the background of the present work with 
some short reviews on related works, basic concepts and prepare basic tools needed in the analysis of the rigidity. 
In sec. mil we develop our scheme to compute the rigidity of structural glasses and apply it to the case of a binary 
mixture of soft-spheres. Finally in sec. IIVI and FVl we summarize the preset work and discuss some related problems. 
In the appendices we present some technical details. The brief account of the present work has been reported in [49|. 

II. THE BACKGROUND 

In this section first we discuss some basic aspects of response to externally induced shear deformations in sec. Ill Al 
Then in sec. IIIBI we discuss a generic mean-field picture which suggests an intimate connection between the order 
parameter and rigidity in solids. Finally we review the cloned liquid method in sec. Ill CI 

A. Linear response to shear 

1. Response to shear : a paradox and a lesson 

A material is said to be in a solid state if its rigidity (shear-modulus) /i is positive. Let us consider N particles put 
in rectangular container which can be deformed by simple shear. Under a simple shear of shear-strain 7, the volume 
V (and thus the number density p ~ N/V of the particles) remains unchanged and only the shape of the container 
changes. (See. Fig.[T]) 



Let us consider how much change occurs in the free-energy F of the system by the simple shear. Assuming that 
the free-energy F = N f is an analytic function of the shear-strain 7, we can expand the free-energy par particle 

/(7) = Fij)/N as, 



/(7) = /(O) +<71+\^i + 



where we find the stress a and the rigidity fi as, 
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However, there is a plain fact that thermodynamic free-energy (par particle) / should not depend on the shape of 
the container (remember that the density remains unchanged under the simple shear deformation). 
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lim F{'y)/N = const. 



(3) 



Thus we conclude that rigidity defined in the thermodynamic sense must always be zero, whether the system is a liquid 
or a solid j50| ! Apparently this goes against to our basic intuition that in solids /i > 0. 

The paradox described above suggests that the thermodynamic limit N ^ 00 and small shear-strain limit 7 — > 
do not commute in solids: linear response theory which is built in the limit limjv->oo linr-y-^o fails to predict what 
actually happens in the thermodynamics lim^_).o limTv-s-oo ■ The breakdown of the commutation of the two limits is a 
consequence of the breaking of the transnational symmetry in solids. Physically the break down of the linear response 
theory means that as a system becomes a solid, not only the elasticity but plasticity must emerge simultaneously. 
The plasticity, which means non-linear responses like yielding or stress drops, recovers the translational symmetry at 
the macroscopic level required by the thermodynamics. In this respect, the idealized elastic bodies which appear in 
macroscopic continuum descriptions are purely hypothetical and thermodynamically unsound objects. In the present 
paper we actually limit ourselves to the linear response theory but we shall always keep in our mind the fact that the 
linear response theory must fail. 

Finally let us note that we strictly consider shear-strain control protocols instead of shear stress control protocols in 
the present paper. In the latter case, static formulations are impossible and one essentially studies the rheology where 
one finds that "everything flows under shear" . Physically the fact that even perfect crystals flow [5l[ is intimately 
related to the static paradox discussed above and the plasticity is the mechanism of the flow of solids under shear. 



2. Static fluctuation formula 



Let us consider a generic system of N particles (i = 1,2,...,A^) of mass m at position r.^ = {xi,yi,Zi) in the 
laboratory frame. More generally we denote components of a position vector r^ in a d-dimensional space as x^ with 
/i = 1, 2, . . . , d. We assume the Hamiltonian is given by. 
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i<j 



where r^- = r^ — r^ and 



\vij\ is the distance between the i-th and j'-th particles. The term K is the kinetic 



energy with pi being the momentum of the i-th particle. For the potential energy U we assume the simple two body 
interactions just for the sake of simplicity. 
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FIG. 1. Schematic picture of the system under static shear, a) The mobile particles (open circles) are bounded by "wall 
particles" (filled circle), b) Laboratory frame {x,y,z) and sheared frame {x',y',z'). 



The free-energy of the system at temperature T can be written in terms of the partition function Z of the eanonieal 
ensemble as, 



pF(T,V)^ log Z 



1 
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(5) 



t<j 



where /3 = l/k^T is the inverse temperature, A is the thermal de Brogile length A = h/ \/2'Kmk^T with fee being the 
Boltzmann's constant. In Eq. ([5]), the symbol V is meant to represent symbolically the integration volume, i. e. the 
container of the system including not only its volume V but also its shape, which will be important in the following. 

In order to study the rigidity against simple shear deformation, we consider a container with two boundary walls 
which are normal to the z-axis and separated from each other by distance Lz as shown in Fig[TJ To impose a shear- 
strain 7 on the system, we simply displace the top wall by an amount 7L2 into x-direction. The shear-deformation 
changes the boundary of the integration volume in Eq. ([5]). Let us denote the deformed boundary symbolically as 
V(7). However, note that the volume of the system V and thus the number density p ~ N/V remain constant under 
this simple shear deformation. 

Now we wish to find an expansion of the free-energy in power series of the shear-strain 7 as Eq. ([T]) . To this end, it 
is convenient to change the integration variables to those of a sheared frame with x' ,y' and z' which are related to the 
original laboratory frame as, {x, y, z) = {x' -\-z'^, y', z') (See Fig.[T]). Note that the boundary of the original integration 
variables V(7) evolves with the shear-strain 7 but that for the new variables remains the same as the original one 
V(0) = V. But now we have to express the Hamiltonian, which is written in terms of the original coordinates, by the 
new ones. This can be done by simply expanding the Hamiltonian in power series of the shear-strain 7 assuming that 
it is small. In practice it is useful to notice. 
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where xij = Xi — Xj and Zij = Zi — Zj. As the result, one easily finds a microscopic expression of the shear-stress 
defined in Eq. ©, 
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where Xij = {xi — Xj)/rij and i^ = (zi — Zj)/rij, and v^"' ^ d"v{r)/dr^ . Similarly the explicit expression of the 
rigidity defined in Eq. ([2]) is found as |26| . 
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In Eq. ([5]) (...) denotes the thermal average evaluated with zero shear-strain 7 = 0, 
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(10) 



The formula Eq. ([8]) is the static fluctuation formula of the rigidity [26|. The 1st term b on the r. h. s of Eq. (|8]) 
is the so called Born term |3a |. It represents the affine response of the system against shear, which is finite even in 
simple liquids. The 2nd term is the non-ajfine correction term due to stress relaxations which play crucial roles as we 
discuss in the present paper. Let us call it as the fiuctuation term in the following. The difference of the nature of the 
affine and non-affine terms will become more clear shortly by discussing the dynamic response against shear. It may 
be instructive to note that quite analogous fluctuation formula is known for the rigidity of spin-waves in ferromagnets 
or the helicity modulus of superconductors j52| where the non-affine part is represented by correlation functions of 
spin or super currents. 

In liquids, the rigidity /i vanishes. It means an exact cancellation of the Born term and the fiuctuation terms. The 
exact cancellation may look somewhat surprising but it is a direct consequence of the translational symmetry in liquids. 



The vanishing of the rigidity in hquids can also be seen formally by considering Mayer expansions, i. e. expansions 
in power series of the number density p = N/V, of the free-energy of liquids. Such an expansion should converge in 
liquids just like high temperature expansions are convergent in paramagnets. At each order in the expansion, one 
finds finite sized clusters. An infinitesimal shear-deformation, which is just an infinitesimal change on the boundary 
condition, should not change the contribution from the finite sized cluster to the free-energy per unit volume in the 
thermodynamic limit V^ — > oo. Thus the rigidity is zero at each order of the density in the thermodynamic limit. The 
vanishing of the rigidity at the level of linear response means the two limits y — >■ cxj and 7 — >■ rfo commute in liquids 
(see sec. Ill A 1[) . Similarly one can consider high temperature expansions in spin systems by which one arrives at the 
same conclusion that the spin-wave stiffness must be zero in paramagnets. 

Finally let us note that it is possible to obtain the zero temperature limit of the fluctuation formula of the rigidity 
Eq. ^ within a harmonic approximation as we discuss in sec. Ill A 71 It is very instructive to note that the vanishment 
of the rigidity in the unjamming transition of a class of systems with contact forces at zero temperature is due to an 
exact cancellation of the Born term and the non-affine correction term WM ^3] . 



3. Dynamic fluctuation formula 

Now let us briefly discuss more general, dynamic fluctuation formula for the linear response to shear-strains which 
describes the elasticity and viscosity in a unified manner. Although our focus is put on the elasticity rather than the 
viscosity in the presen t paper , the overview will become useful in discussions. Such a unified view point has been 
developed also in Ref. |5J-l5g . 

Linear response against shear deformations can be seen experimentally by measurements of the linear- viscosity (See 
for example |57H6(j |). In sec. Ill All we emphasized that for arbitrary small but finite strength of perturbation 7, 
static linear response break down in the N ^ 00 limit in solids. Similarly dynamic linear response theory will break 
down at large enough time scales by overwhelming non-linear contributions. However, it is still possible to delay the 
dominance of the non-linear responses and bring it out of a given observation time (or frequency) window by choosing 
sufficiently small shear-strain 7 |58| . Of course one has to choose smaller 7 for larger time window. 

Within the equilibrium linear response theory the shear-stress a{t) at time t can be related to small changes of the 
shear-strain 7(i') in the past t>t'> —00 as (54| . 

6a{t) = bj{t) (3 f^ dt' ^^f/\ {t'). (11) 

Here b is the born term defined in Eq. ^. The correlation function Ca{t,t') is the shear-stress auto-correlation 
function defined as, 

C^{t,t') = N{a{t)a{t')) (12) 

We show in Appendix \K\ a simple derivation of Eq. ((TT|) . The 1st and 2nd term on the r.h.s. of Eq. pT|) represents re- 
spectively the instantaneous rigidity, which describes the affinc part of response to shear, and the non-affine correction 
term due to shear-stress relaxations. 

By doing an integration by parts in the r. h. s. of Eq. (|lip we can find an alternative expression for the linear 
response, 

Sa{t) = {b- N/3{a^))-fit) + Np{af -/{~oo) +p j dt'C^it, t')j{t') (13) 

where 7(t) = d'-f{t')/dt' is the shear-strain rate. 

More specifically, let us consider two typical experimental protocols, 

• Step like shear : 7(t) = 7^(i) 
The stress relaxation after then step like perturbation can be found to be, 

^(i) = <5(T(t)/7 = 5-/?(a(t,i) - a(<,0)). (14) 

The initial value is nothing but the Born term 

/i(0) = 6, (15) 



and it relaxes down to the static rigidity /i defined in Eq. ([8]) in the large time limit, 

^ = lim n{t) = b - NI3{{a^) - {af). (16) 

Thus the remnant stress under external strain beeomes finite in solids and vanishes in liquids at the level of the 
linear response. 

Furthermore, in the case of liquids the rigidity vanishes /i = 0. This implies (see Eq. (fT6|)). 

h = Nli{(j'^). (17) 

Here we used (<t) = which must always hold in liquids. It is worth to mention that the relation Eq. ()17p 
means the born term (instantaneous rigidity) and variance of the thermal fluctuation of the stress are intimately 
related in liquids in equilibrium. As a result, the stress relaxation Eq. (J14p becomes directly proportional to the 
shear-stress auto-correlation function, 

M(i)=/3C.(t,0). (18) 

• Constant shear-rate : 7(t) = 7 
Let us consider here liquids for which the rigidity is zero ^ ~ 0. Then from Eq. (|13p we find, 

/•oo 

5(7/7 = 77 r^ = P dTG,{T) (19) 

Jo 

where rj is the shear-viscosity. The last equation is nothing but the Green-Kubo formula for the shear- viscosity 

Usually we tend to think that rigidity and viscosity are very different properties. In this respect it is very interesting 
to note here that both the non-affine correction term of the rigidity, which appears in the 2nd term on the r. h. s of 
Eq. p6)) . and the viscosity given by Eq. (fT9|) are related to the spontaneous thermal fluctuation of the shear stress. 
Quite remarkably recent numerical studies of glasses have shown that the non-affine correction term is significant even 
in the zero temperature limit [62, Ifij] . (The zero temperature limit of the fluctuation formula of the rigidity Eq. ^ 
is given by Eq. p7p ). This implies rigidity and viscosity are related in glasses. 

4- Visco-elasticity 

Supercooled liquids are visco-elastic materials because of the two-step relaxation: the shear-stress auto-correlation 
function Ca-ir) or equivalently the stress-relaxation ^(t) (see Eq. ([T5|) ) exhibit plateau like behaviours as shown 
schematically in Fig. [2l 

The stress relaxation function fi{T) starts from the instantaneous rigidity, i. e. the born term b at time r = (see 
Eq. (|15p ). Very at the beginning it exhibits a rapid decay as indicated in Fig. [51 which we call as 'initial relaxation'. 
The duration of the initial relaxation is expected to be of the order of the vibration time tq = ^mo? jt with a being 
the microscopic unit length scale and e being the uni t energ y (or collision time p^^ I "^ ^frajk^ ). The rapid initial 
relaxation can be easily seen in numerical simulations [68 70] . Whether it lies within the time window of experiments 
would depend on the specific systems under study [57lh60i] . 

Then the relaxation slows down significantly exhibiting a plateau like behavior [60, l68i - l7C)| . which we call as '/3- 
relaxation'. Let us call the height of the onset of the plateau like region as Goo in the present paper. If the plateau is 
ideally extended for a long enough time scale, the system would appear as a solid with a finite rigidity (see Eq. ([T6|)), 
which is the main concern of this paper. 

The temporal solid state is characterized by the cage structure: the particles are collectively confined by the 
surrounding particles within a narrow space called as cage. Within the framework of the RFOT, including both the 
MCT [ll|, [i3l and cloned liquid approaches [l^ [20| , the size of the cage A is viewed as the order parameter which 
indicates the presence of amorphous solid state (see sec. IIIC ip . 

Finally the a relaxation starts at longer time scales due to reorganization of the cages by which the residual stress 
^(r) eventually decays down to so that the viscosity 77 Eq. p^ can be defined. As the temperature is lowered the 
relaxation time of the supercooled liquid Tq, increases rapidly and eventually exceeds the experimental time scale. As 
the result the system falls out of equilibrium in experimental time scales. 

Still visco-elastic properties at lower temperatures can be studied systematically by observing aging effects [SJ, UM 
which allow us to explore the glassy phase below Tk. For example, a typical stress-relaxation protocol [53 . [60| goes 
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FIG. 2. Schematic, idealized picture of stress relaxation of supercooled liquid/glass under aging (left) and an intuitive free- 
energy landscape picture (right). The stress relaxation function fj,{T;tw) is presented in a unit such that /i(0) = 6=1. Here 
b is the Born term [33] which represents the instantaneous rigidity of system. The time r is presented in a unit such that the 
microscopic time scale tq, corresponding to typical vibration time of the particles, becomes 1. The waiting time fw increases 
from the left to the right curves. If the relaxation time t^ is finite, the curves converge to a limiting (equilibrium) curve at 
large enough waiting time iw 3> t^. We consider that the stress relaxation consists of three regimes. Within the initial very 
short time regime whose duration is just of the order of the unit time scale tq, the stress relaxes rapidly starting from the 
instantaneous rigidity b down to some value, which is often called (somewhat misleadingly) Goo- Most likely this is due to 
fast modes of purely elastic relax;ation in the vicinity of inherent structures (IS) [6J, [Sg and that Goo roughly corresponds 
to rigidity of inherent structures /xis. In the subsequent /3-relaxation regime, the stress value remains within a plateau-like 
region in the ran ge G p < ^{t) < Goo- We argue that this relaxation takes place within metastable states, each of which is a 
metabasin (MB) [6y, [631 i- s. a union of ISs with large enough mutual overlaps rather than individual ISs. Furthermore we 
argue that Gp corresponds to rigidity of metastable states fi which we analyze in detail in the present paper. The reduction 
of jl with respect to the rigidity of inherent structures pis is attributed to thermal excitations within metastable states, such 
as mutually independent, localized plastic deformations. Thus the gap between Gp and Goo is expected to become smaller 
at lower temperatures. Finally the stress relaxes down to by the «-relaxation by which the system explores larger phase 
space. While the a relaxation strongly depends on the waiting times fw, the initial and /3 relaxations depend on iw only weakly 
(quasi-equilibrium) . 



as follows: 1) the system is prepared in the equilibrium liquid state at a high enough temperature above the glass 
transition temperature Tg , 2) the temperature of the bath is quenched to a target temperature T below Tg, 3) the 
system is left in a quiescent state for a waiting time t^ after which 4) a finite shear-strain of amplitude 7 is applied. 



s-fit) = je{t-u). 



(20) 



Again the stress relaxation starts from the instantaneous rigidity, the born term b and exhibits the two step relaxations 
(-1- initial relaxation) as shown in Fig. [21 But now the stress relaxation depends not only on the time r elapsed after 
switching on the shear-strain but also on the waiting time iw thus we denote it as /L((T;iw)- Most notable feature is 
that a-relaxation becomes strongly dependent of the waiting time iw while the /3-relaxation only weakly depends on 
the waiting time fw. Then it would be meaningful to define a limiting relaxation curve for the /3- regime. 



Mcq(T) = lim M(T;iw)- 



(21) 



As suggested by the studies of aging effects, in particular the dynamical mean-field theory of spin-glasses [74I, it 
is natural to expect that the fluctuation formula Eq. p^ still holds in the /3-regime and that a slightly modified 
fluctuation formula applies for the a-regime (See Appendix [A| . 

The system will behave as a solid with finite rigidity as long as the time scale lies within the plateau-like region. 
However there is an ambiguity to define such a rigidity because the plateau is not completely flat at finite time scales. 
As noted above we call the value at the onset of the plateau-like region as Goo in the present paper. Physically it is 



more interesting to examine the value of the plateau in the asymptotic limit (Sec Fig. [2]), 

Gp = lim ^eq(T) (22) 

Note that as long as the system is ergodic, which is certainly the case at high enough temperatures in the liquid phase, 
the limiting curve Eq. (pij) covers not only the /3-relaxation but also the a-relaxation by which the shear stress relaxes 
down to 0. Thus for this Gp to become non-zero, there must be ergodicity breaking, i. e. the ideal glass transition. If 
the putative ergodicity breaking take place at some temperature, say Tne, then the asymptotic value of the plateau 
at Tne can also be obtained following the equilibrium liquid state at T > Tne as, 

Gp(Tne) = Um lim Aieq(T) (23) 

It is not known yet whether such a true ergodicity breaking (ideal glass transition) really exist in real systems. 
However even if it does not exist in the rigorous sense, we may take a pragmatic view and look for some crossover 
height Gp between more j3-like regime and more a-like regime. 

5. Free-energy landscape picture 

Let us discuss the visco-elasticity of supercooled liquids within the phenomenological free-energy landscape picture 
shown schematically in Fig. [2] [23, |64| - [67| . The following arguments are based on simple intuitions on some basic 
energetics so that they are admittedly naive in many respects and thus they should be taken with a grain of salt. Here 
we focus on deeply supercooled, low temperature regime. In sec. Ill A 91 we briefly discuss the free-energy landscape 
at higher temperatures. 

The most basic ingredients in the free-energy landscape at low temperatures are local energy minima which are 
called as inherent structures (IS) [6J,|65|- In the phase space, each inherent structure is surrounded by a local region 
called as basin whose energy landscape can be described well by a purely elastic Hamiltonian (See Eg ([Ml) below) . 

A set of inherent structures and their basins may be grouped together as a metabasin (MB) [60, [63] • From the 
point of view of the mean-field theory which we discuss in sec. Ill CI a metabasin rather than an inherent structure 
should be identified as a metastable state which are pure states in the mean-field theory [73| . Within the mean-field 
theory, metastable states are characterized by the cage size A (See Eq. (|6T|) ). Two difi'erent inherent structures, say 
ISi and IS2 are associated with a common metastable state if the average mean-squared distance is smaller than the 
cage size, 



2Nd 



J2(r?'-r?n'<A (24) 



Here r^ ^ and r,- ^ are configurations of the particles at the two ISs. For example consider two sets of particle 
configurations such that one configuration is created from the other by a localized plastic deformation involving 
just a few particles. Differences between the two configurations on the rest of the system will be just smooth 
relative elastic deformations whose amplitude decay to at long distances from the region where the localized plastic 
deformation took place [23|. Such a pair of configurations, both of which are energy minima, can be easily generated 
explicitly for example by athermal quasi-static shear (AQS) processes [33. \33i. I&a. [74l - l77| . In the thermodynamic limit 
N —)' 00 the mean-squared distance between such a pair of states which appears on the left hand side of Eq. (|M)) 
will vanish. Although they are distinct inherent structures which are not accessible to each other by continuous 
elastic deformations, so that each one is out of the basin of the other, they should be regarded as belonging to the 
same metastable state since their thermodynamic properties must be the same. This is an extreme example but the 
argument can be repeated for the cases involving two, three, ... localized plastic deformations which are independent 
from each other. 

Within a given metastable state or metabasin, there must be an inherent structure which has the lowest energy, 
i. e. the ground state within the metabasin. All other inherent structures are excited states in the metabasin. As 
the temperature increases, not only the ground state but also the excited states become more populated. Moreover 
different metastable states may merge together by raising the temperature. As the result the entropy of the metabsin 
increases with the temperature. 

We consider that the stress relaxation of deeply supercooled liquids may be described within the free-energy land- 
scape picture as follows (see Fig. [5]). First, we consider that the initial rapid relaxation is mainly due to the dynamics 
in the basins of individual ISs, which arc dominated by fast modes of the purely elastic dynamics which we discuss 
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in sec. Ill A 71 The subsequent slower process in the plateau-region, i.e. the /3-relaxation would involve possibly two 
different mechanisms. One is due to slower modes of the elastic dynamics around inherent structures. The other 
is due to thermally activated transitions between different inherent structures in a common metabasin (metastable 
state). The latter is necessarily inelastic (plastic) deformations as we discuss in sec. Ill A 81 Finally the a-relaxation 
may be attributed to thermally activated transitions between different metabasins (metastable states). 
The above phenomcnological picture implies thermal averaging may be factorized as, 

(...) = [(... )ms1 (...)ms^[(--->is]mb (25) 

Here (. . .)ms ^^nd |. . .] represent respectively thermal averaging within a metastable state and over different metastable 
states. Similarly (. . .)is and [. . .]mb represent respectively thermal averaging in the vicinity of an inherent structure 
(IS) and over different ISs within a common metabasin (MB) (or metastable state). Then thermal fluctuations of a 
physical quantity, say O, could be decomposed into intra-state and inter-state fluctuations as, 

iO') ^ {Of = l{0\,s] - l{0)usf = 1{0')ms - {Ofusl + l(0>^s] - I(0>MsF • (26) 

■^ . ■' , ^ 

intra-statc fluctuation (initial+/3) inter-state fluctuation (a) 

Quite remarkably the cloned liquid method enables disentanglement of the intra-statc and inter-state fluctuations and 
decompose thermal fluctuations into the two parts as Eq. ((26)) as we discuss later in sec III C 41 

Moreover, the intra-state fluctuation itself may be decomposed into thermal fluctuations around inherent structures 
and those between different inherent structures in a common metabasin, 

(02)ms - {0)ls = [{0')is]mb - [(0)is]^B = [(O')is - {0)!s]mb + [{0)Umb - [{0)is]Iab (27) 

-^ , ^ ^ , ^ 

fluctuation around IS (initial) inter-IS fluctuation (/3) 

Let us note that /3-relaxation may not be entirely due to the thermal activation between different inherent structures 
but partly due to slower modes of the elastic fluctuations around inherent structures as discussed above. How the 
two different relaxations mechanism are mixed in the dynamics is not clear at the moment. 

6. Hierarchy of rigidities 

The hierarchy of the thermal fluctuations naturally implies hierarchy of rigidities. The rigidity as defined by the 
static fluctuation formula Eq. ([8]) (or Eq. (|16|))) contains non-afhne corrections represented by the correlation function 
of the shear-stress fluctuations. By decomposing the latter into several pieces following the above prescription we 
find, 

fi^lib- N/3 ((a2)is - {a)fs)]MB -N(3 {[{a)fs]MB - [{<j)is]lm)]] -NP ([(a)Ll - I^Msf ) (28) 

Mis 

A 
where we introduced the rigidity of inherent structures /ijs , 

^ilS = h-Np{{a\s-{cy)h) (29) 

and the rigidity of metabasin or intra-state rigidity ft, 

A = 5-iV/3((a2)Ms-(^)Ms) (30) 

with b being the Born term Eq. ([9]). 

We speculate that the rigidity of inherent structures /iis and that of metastable states (metabasins) /I corresponds 
respectively to the upper and lower bound of the plateau region (/3-relaxation) of the shear-stress relaxation (see 
Fig.©, 

Goo ^ m = Np{a)\s (31) 

Gp ^ /i = Np{a)l,^ (32) 

Since the initial and /3-relaxations weakly depend on whether the system is in equilibrium or out-of equilibrium, we 
expect that the quasi-static evaluations of Goo and Gp are valid for the supercooled liquid in equilibrium as well as 
glasses out-of equilibrium. 
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The last equations of Eq. ((3T|) and Eq. (|32| relates the rigidity at different levels to the variance of the shear-stress 
fluctuations at the corresponding levels. It is easy to see that they must hold in liquids. First note that the total 
rigidity must vanish /^ = in liquids. Then note that the stress must vanish after fully taking the thermal average 
("') — [(o')msI — 0- Using these two in the expression Eq. ([^5)) one immediately finds the last equations of Eq. ([31]) 
and Eq. ([32|) . Note also an analogous relation Eq. (|17p which reads 

b^ H{t = 0) = NP{a^). (33) 

It relates the born term b (instantaneous rigidity) to the variance of the total thermal fluctuation of the stress {<j'^). 

7. Rigidity of inherent structures 

Here let us discuss more explicitly the rigidity //is of inherent structures (IS). An IS [6J| is a configuration of the 
particles in an energy minimum around which the particles make purely elastic fluctuations. In the close vicinity 
of an IS with the particle configuration {(ris)i} {i = 1, 2, ... , N), the potential part of the Hamiltonian Eq. ((4]) is 
approximated as, 

U = U{{{x'^)is}) + Uci{{{u'^)is}) + unharmonic corrections. (34) 

Here C/ei({(uf )is}) is the harmonic part of the energy, 

[/ei({«)is}) ^Ij2Y. ^r<".' < = ^'' - (^is)r (35) 



with the matrix H being the Hessian matrix, 

(36) 



^r = '^M-'^ J2 ^'"'in,k) - v^'^in,) v^'^ir,,) ^'^("'^") 



dx^dxj 



IS 



where .-.Ijg means to take derivatives at the IS, namely x^ = (a;is)f. Then the intra-IS rigidity Eq. (|29p can be 
expressed fully microscopically [6J, uE, UM ^^i 

1 ^ 
^^ls=b--Y.^^■{H-'m^, (37) 

where H~^ is the inverse of the Hessian matrix Eq. ([55)1 and the H is defined as, 

H. = V,a|is (38) 

with a being the shear-stress Eq. ([7]) and V; = {d/dxi,d/dyi,d/dzi, . . .). More precisely, the born term Eq. ^ at 
the IS can be written as, 



^=^E^5'(^is).(^is)„ (39) 

and the H at the IS can be written as. 



TV 



Sf^E^rC^is)/. it^-x,y,z). (40) 

k 

Within the harmonic approximation to the potential Eq. p4p , the fluctuation formula of the rigidity at finite tem- 
peratures Eq. (|5]) reduces exactly to Eq. (|57)) [78| . Non-harmonic corrections which can be treated perturbatively 
amount at most to 0{T) terms which vanish in the limit T — )• 0. 

It is instructive to rewrite the fluctuation formula Eq. (|37p in terms of the cigen modes of the Hessian matrix as 

"2 



^--^-^E%- (41) 



N '^ X 
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Here A's are the eigcn values of the Hessian matrix which satisfy the equations, 

E^rf^O-.-),A = ^%/^),A (42) 

0» 

with U(j^^^)\ being the normahzed eigen vectors. The eigen values A of the Hessian matrix at the inherent structures 
(energy minima) are not negative by definition. In Eq. (|4ip we also introduced the projection of the S field onto the 
eigen modes, 

SA = E(t/t),,(,,)Sr. (43) 

The shear-stress relaxation ^{t) can be studied easily within the harmonic approximation. For simplicity we assume 
the equations of motion with viscous damping, 

^^^^^^ = -T.H^,v^ (44) 

with 771 and 77 being the mass of a particle and shear- viscosity. The coupled equations of motions can be transformed 
to independent equations of motions for the eigen modes u\ = ^j-,- u)(t^^)A,{i.M)^i'' 

where we introduced characteristic frequencies 10 defined as. 



L0{\) = y/Xj^l (46) 

We further simplify the problem considering the over-damped limit neglecting the inertial term. Then the shear-stress 
relaxation jiit) after the step like perturbation is obtained as, 

u(t) = b- dwD(w)^^[l -e-("'/")*]. (47) 

Jo muj^ 

where we introduced the distribution function D{uj) of the characteristic frequencies w, 

^H^^E'^(^-v/V^)- (48) 

A 

Note that in the large time limit t — ?► 00, fi{t) converges to the static rigidity //is given by Eq. (|4ip as it should. 

The spectrum of eigen values of the Hessian matrix of the ISs is known to constitute a continuous band extending 
down to w = so that the stress relaxation fi{t) cannot be a simple exponential relaxation even in the harmonic 
approximation. However convergence of fi{t) to the asymptotic value /iis would be not too slow if the integrand of 
Eq. (j47p is not singular in the w — > limit. Indeed the numerical analysis of inherent structures of a model structural 
glass system in [T^ (See Fig. 5 of the latter reference) suggests this is the case such that the a large fraction of the 
whole reduction of the rigidity from b to /iis would be attained within the order of the microscopic time scales. This 
observation is consistent with our view that the initial rapid decay of stress-relaxation (see Fig. [5]) is most likely purely 
elastic dynamics in the close vicinity of ISs. 

8. Rigidity of metabasin 

In sharp contrast to the purely elastic deformations discussed above, the transitions between different ISs, even 
within a common MB (See Fig. [5]), necessarily involve some inelastic {plastic) deformations j27i] . It is natural to expect 
that they are thermally excited plastic deformations, which are mostly independent from each other. We consider 
such processes, together with slower modes of the elastic dynamics, contribute to the slower relaxation around the 
plateau, i. e. the /3-relaxation |65l |. 

Let us come back to the example of the pair of IS which we discussed at the beginning of sec. Ill A 51 They are very 
closely related to each other such that the transition between the two involves only a localized plastic deformation 
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^^(7)) 



FIG. 3. Schematic free-energy landscape of metabasin under external strain and corresponding static stress-strain curve. The 
free-energy f (7) (bold line) and stress {cr(7)) = dF{'y)/d'y (bold line ) is plotted against the 7-axis. In this example two 
inherent structures (IS) with rigidity /^is, represented by the parabolas, belong to the metabasin. At zero temperature, the 
level crossing between the two result in a discontinuous stress drop in the stress strain curve. At finite temperatures, the 
singularity at the level-crossing is removed such that the renormalized rigidity of the metabasin ji can be defined. 

involving a few particles. Such transitions between two ISs are already beyond the scope of purely harmonic theories. 
In such a process, typically an eigen value of the Hessian matrix becomes negative at the transition [6j]. This means 
that elasticity as defined in the most strict sense breaks down immediately moving just a little bit away from the close 
vicinity of an IS. 

However it is reasonable to expect that some sort of renormalized harmonic theory which takes into account the 
effects of some weak plasticity should exist at macroscopic level. In this respect it is instructive to recall, for example, 
the fact that the rigidity of spin waves of ferromagnets can still be defined in the presence of some small densities of 
topological defects such as pairs of vorticies (dislocations). In the latter case, thermally excited topological defects 
renormalizes the the rigidity down to a smaller value as the temperature is increased to higher temperatures and 
eventually enables melting of otherwise elastic system [52|, [T^ • 

Our picture for the rigidity of metabasin can be summarized as the following. Different ISs within a common 
metastable state or MB have large enough mutual overlaps in the sense of Eq. (P^ . They are likely to be related to 
each other by localized plastic deformations [23, [S^, S, |6J, (TJ, [T^, ISOl ■ Among all ISs within a common MB there is 
an IS which has the lowest energy, i. .e the ground state within the MB. All other ISs in the same MB are excited 
states, which can acquire larger statistical weights at higher temperatures. In a sense the dynamics within a MB is 
"reversible" [23, [Sa, [631 because a transition to an excited state would be often followed by a return to the ground 
state. From Eq. (P^ and Eq. ([50]) we find the intra-state rigidity jj. is reduced with respect to the rigidity of IS /iis 
due to the stress fluctuations between different ISs in the common MB, 



M 



= m - NI3 ((((7)?s>MB - ((a)is>U) • 



(49) 



In our interpretation, the non-affine correction which appears on the r. h. s. is plastic in sharp contrast to the elastic 
non-affine correction in the basins of inherent structures discussed in sec. Ill A 71 In the low temperature limit, 
however, the thermal excitations will be suppressed so that it is natural to expect that the two rigidities become the 
same. 



lim u 



lim ^is- 



(50) 



In order to illustrate our picture in an alternative way, we show an intuitive free-energy landscape picture of a 
metabasin plotted along the axis of external strain 7 in Fig. [3] together with the corresponding static stress-strain 
curve. Local energy landscape of the basins around each inherent structure should be described well as a parabola 
with curvature /Uis- The positions of the bottoms of the parabolas of different ISs would be slightly displaced from 
each other. At a given 7, the ISs can ranked according to their energy: ground state, 1st excited state, e.t.c. Along 
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the 7-axis there can be level crossings among the ISs, i. e. plastic deformations. At zero temperature each level 
crossing would be reflected on the stress-strain curve as a stress-drop. However at finite temperatures, the /ree-energy 
landscape and the corresponding stress-strain curve should be rounded. Then it is natural to expect that the /ree- 
energy landscape of the metabasin as a whole may be described well as a large parabola with curvature /t which is 
reduced with respect to nis at finite temperatures. 



9. Free-energy landscape at higher temperature 

Let us finally turn to discuss on the free-energy landscape at higher temperatures which will provide us useful 
insights on how glasses melt. Within the RFOT [5|-Q, the critical behaviours predicted by the mode coupling theory 
fMCT)[lll. [l2| are considered as reflections of melting or emergence of amorphous metastable solid states. If this is 
the case there must be an intimate analogy between the MCT criticality and the spinodal criticality of superheated 
crystals. In sec. Ill B I we will discuss implication of this point of view on the rigidity problem. 

It is important to note however that, at variance to the spinodal criticality of superheated crystals, the MCT 
criticality is relevant not only in the low temperature side but also in the high temperature side above T^. Actually 
MCT itself found the criticality in the higher temperature side. This point is very important in supercooled liquids 
which exhibit the two-step relaxations as the one shown in Fig. [5] not only at temperatures lower than Tc but also at 
higher temperatures around T^. 

The key notion to describe the free-energy landscape of supercooled liquids at T — > T+ , corresponding the inherent 
structures for the lower temperature side, are saddles. This has been substantiated by extensive studies of a class of 
solvable mean-field spin-glass models and realistic model glass forming systems [^, |4j, |81| . The basic result goes as 
follows. The dynamics of a super-cooled liquid at T — >■ T+ can be viewed as a point in the phase space going through 
a narrow channel along which the energy can be decreased. As the time goes on, the number of eigen modes of the 
temporal Hessian matrix with negative eigen values becomes less and less. In the large time limit the system ends 
up in a marginally stable state: all the eigen values of the temporal Hessian matrix turns out to be positive but the 
lowest one is touching zero. According to this picture, the main mechanism of the slow dynamics as T -> T^ is a sort 
of critical dynamics distinct from the activated dynamics which becomes relevant at lower temperatures. 

The essential difference between saddles and inherent structures are due to the slow modes with smaller eigen 
values. The saddles which are relevant for the long time dynamics as T ^ T+ contain a large number of positive 
eigen values and a small number of slightly negative eigen values. Then it is tempting to speculate that a somewhat 
similar elastic description of the stress relaxation as expressed by Eq. (|T7)) would be also possible for T > T^ for short 
enough time scales. 

Conversely if a supercooled liquid initially at a temperature below Tc is heated up to a temperature above T^, the 
temporal Hessian matrix will acquire a fraction of negative eigen-values. Even with one negative eigen value, the 
system will yield by arbitrary weak shear perturbations leading to collapse of a metastable amorphous solid state. 
This suggests that the static rigidity of a (temporal) metastable amorphous state would vanish approaching T^ from 
below. 



B. Connection between the order parameter and rigidity : mean- field points of view^ 

As we discussed in sec III All emergence of rigidity /i > in solids means break down of the commutation of the 
two limits, namely the thermodynamic A^ — > oo and small shear-strain limit 7^0. Thus the rigidity detects the 
liquid-solid phase transition as sensitively as the order parameters. Here we show that several mean-field theories 
generically suggest an explicit connection between the rigidity and the order parameter. 

In the introduction of the present paper we mentioned the Horn's rigidity crisis scenario |37| which says that 
melting of solids may be signaled by vanishing of the rigidity /i. The problem is how it vanishes: continuously or 
discontinuously? Usually Horn's conjecture is interpreted as suggesting continuous vanishment. This is an interesting 
and much debated question in the case of superheated crystals |38l - l4l| approaching their metastability limit at Tg 
[82| . The same question is relevant also for supercooled liquids and glasses since, as emphasized in a recent work j42| , 
the dynamical transition at T^ (or the MCT critical temperature T^) can be regarded as an analogue of the spinodal 
behaviour of superheated crystals at Tg. We show below that some mean-field theories indeed suggest a common, 
universal behaviour of the rigidity approaching the metastability limits of solids including both the crystalline and 
amorphous systems. 

In the following we first develop a simple mean-field description of the rigidity of a superheated crystalline system 
(ferromagnct). Then we recall the predictions by the MCT on the fate of the rigidity of amorphous solids. 
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FIG. 4. A model Landau potential. The potential is defined in Eq. (|5ip . The temperatures are T = 5/8, 6/8, . . . , 10/8 from 
the bottom to the top curves. The curves at the spinodal temperature Tg = 1 and the 1st order phase transition Tp = 3/4 are 
indicated. The two minima on the left and right corresponds respectively to the paramagnetic and ferromagnetic state. The 
equivalent effective potential of structural glasses and related spin-glass systems is the so called Franz-Parisi potential [SJ] for 
which the horizontal axis should be interpreted as the axis for the Edwards-Anderson order parameter qea or the cage size 
A. In the latter case the two minima on the left and right corresponds respectively to the paramagnetic state q = (liquid 
A — oo) and spin-glass state g > (glass A < oo). The two temperatures Tg and Tp corresponds respectively to the dynamical 
transition temperature T^ (or the MCT transition temperature Tc) and the Kauzmann transition temperature Tk. 



1. Rigidity in a GL theory for superheating 



Let us first discuss the fate of the rigidity of superheated crystalline systems approaching their metastability limits. 
For simplicity we consider Ginzburg-Landau (GL) theory for a vectorial spin system which can realize superheated, 
metastable ferromagnetic phase. By slightly generalizing the usual scalar GL theory for superheating |82l |. we consider 
an 0(2) symmetric (j)^ theory described by a free-energy functional. 



F 



= / d'^r 



|(v0y+5(|0P 



aiy) '^Ty-y'^ + — 



(51) 



where 0(r) = (0i (r), 02 (r)) is a continuous vectorial field in a d-dimensional space. The profile of the Landau potential 
(7(0^) is shown in Fig.Hl 

Let us first examine the equilibrium states. The equilibrium states are uniform in space and parametrized by the 
amplitude of the order parameter M = |0| = yj (j)\ + (j)\ ■ At temperatures below the spinodal temperature Tg = 1, 
there is a metastable ferromagnetic state with non-zero amplitude of the order parameter M > 0. Equilibrium 
1st order phase transition from the paramagnetic state M = to the ferromagnetic state il/ > takes place at a 
temperature Tp = 3/4 lower than the spinodal temperature Tg — 1. 

Upon heating the metastable ferromagnet above Tp up to Tg, the amplitude of the order parameter M exhibits a 
discontinuous jump from a finite value M > (metastable ferromagnetic state) to zero Af = (paramagnetic state ) 
signaled by a preceding square-root singularity. 



M{T) - M{Ts) oc yjT~T. 



(52) 



The mixed character of 1st and 2nd order phase transitions is the generic feature of the spinodal criticality [82l |: 
the order parameter vanishes discontinuously much as usual 1st order equilibrium melting transitions but the jump 
accompanies a preceding square- root behaviour as T — > T^ . 

Now let us consider spatial fluctuation of the order parameter. To this end it is convenient to re-parametrize the 
order parameter using 'angular variable' 9(v) as, 



0(r) = |0|(r)(cos6'(r),sin6'(r)). 
Then the free-energy functional Eq. (j5ip can be rewritten as. 



F 



m 



= /A 



|(V0)2 + £(vh) +g(\$\' 



(53) 



(54) 
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with 

^i = c\$\\ (55) 

Note that the last two terms of the integrand in Eq. ([54|) appear also in the usual scalar GL theory. On the other 
hand, the 1st term of the integrand can be naturally interpreted as the 'spin-wave' term. 

An important observation is that the rigidity /i of spin-waves is simply proportional to the square of the order 
parameter at this mean-field level [52l |. In particular, assuming that the amplitude of the order parameter is uniform 
in space and takes the equilibrium value M{T) we find, 



^l{T) - fi{T,) (X c^/t~T. (56) 

using Eq. ([52]) in Eq. ([55|). Thus at the pure mean- field level the rigidity of the spin- wave n also exhibits a discontinuous 
behavior preceded by a square-root singularity as T — ?► T~ . 

Note that the situation is quite different from the case of the usual second order ferromagnetic transition (e. g. 
(/)^ theory) [52|. In the latter case the amplitude of the order parameter exhibit a similar square-root behaviour 
M(T) ex \/Tc — T approaching the critical temperature Tc from below hut with no discontinuity at Tc- Then using 
the latter in Eq. ([55)1 one finds that the spin- wave stiffness scales as ^i{T) ex c{Tc — T). 

2. Mode coupling theory 

The mode coupling theory (MCT) [111 [l^l predicts that the density auto-correlation function Cpir), such as the 
self-intermediate scattering function, exhibits a two-step relaxation with a plateau similar to the one shown in Fig. [2] 
in supercooled liquids 

and that the ergodocity breaking takes place at the MCT critical temperature T^. the a-relaxation time Tq diverges 
as T ^ r+. Physically this means emergence of amorphous solid states with the cage structure. The amplitude of 
the plateau of the density auto-correlation function Cp{T) is related to the cage size A which is the order parameter 
within the RFOT (see Eq. (|M)) V The MCT predicts that the cage size A exhibits a discontinuous jump from to a 
finite value as T ^ r+ followed by a square-root behaviour for T < T^ 



A(T,) - A{T) ex y^T~T. (57) 

Note that the behaviour is very similar to that of the order parameter of the superheated crystalline system Eg . (|52p . 

The cage size A can also be defined in a static matter (see Eq. (|61]) ). In the cloned liquid approach jl8| - |20| , it is 
obtained by extremizing an effective free-energy which is a function of A (see sec lIIC5"|) . The effective free-energy 
corresponds to the so called Franz-Parisi potential [83[ which has essentially the same qualitative features as the 
Landau potential shown in Fig. |4l The cloned liquid approach [20| predicts that the minimum of the potential 
corresponding to the metastable amorphous solid state A < oo disappears approaching the dynamical transition 
temperature Tj from below, i. e. as T — > T^ . The dynamical transition temperature Tj is considered to corresponds 
to the MCT critical temperature T^. The resultant behaviour of the cage size A is the same as that predicted by 
the MCT Eq. ([57)) . Quantitative comparisons of the results by the dynamic and static approaches are currently 
undertaken intensively |22| . 

Within the convectional MCT, the behavior of the the stress-autocorrrelation function Cair) is obtained by pro- 
jecting it onto the density auto-correlation function Cp{T) [6^[8J]. The former is assumed to be simply proportional 
to the square of the latter Caij) ex Cplr)"^ which follows from the so called factorization approximation [SJ]. As 
the result the MCT predicts a similar plateau behaviour for the stress relaxation function fJ-ir) as for the density 
auto-correlation function Cplr) : the asymptotic plateau value of the the stress-relaxation curve /^(r) = I3C(j{t) (see 
Eq. P^ ) exhibits a discontinuous jump from to a finite value G'p(Tc) as T ^^ r+ followed by a square- root behaviour 
for T < Te, 



Gp(rc)-Gp(T)ex V^c-T. (58) 

Note that this is again very similar to the behaviour of the rigidity of the superheated crystalline system Eq. ([56[). 



3. Discussion 

The discussions presented above suggest common features of the rigidity of crystalline and amorphous solids at 
the mean-field level. First of all, the rigidity does look like the order parameter: it is just the square of the order 
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parameter. As the result, it vanishes discontinuously preceded by a square-root behaviour approaching the spinodal 
like melting temperature from below. Although the preceding square-root singularity somehow signals the melting, 
the discontinuous jump of the rigidity disagrees with the usual interpretations of the Born's original rigidity crisis 
scenario [37l - l40| . 

The phenomenological GL approach like the one presented in sec III B ll is simple and instructive by itself. The local 
mean-field (Landau) potential 5(0^) (see Fig. |4]) has an analogous one for glasses, namely the Franz-Parisi potential 
[83| as we mentioned above. Indeed we demonstrate in Appendix [C] that this approach can be naturally extended 
to the case of vectorial spin-glass system using some replica field theoretic formalism, which implies again the same 
generic behaviour at T^: discontinuous spin- wave stiffness preceded by a square-root singularity. 

However in the main part of the present paper (sec IIIip , we will analyze the static rigidity of glasses starting from 
a microscopic level using the fluctuation formula of the shear-modulus (see sec III A2|) rather than working with the 
phenomenological GL like approach. As the result we will obtain the rigidity of the metastable amorphous solid states, 
i. e. the intra-state rigidity fl Eq. (|30p in power series of the cage size A as Eq. (|156p which reads, 

iiiA) = co-ciA + 0{A^). (59) 

with Co and Ci are positive constants which we evaluated explicitly. Thus we do find an intimate connection between 
the rigidity and the order parameter. Based on the result, we will come back to the question on the fate the rigidity 
in sec. HITEl 

Finally, we should remind ourselves after all that the spinodal like criticality only exists at the pure mean-field level 
and that to what extent the signatures of such a criticality remains in realistic systems is an open question [82| . 

C. Cloned liquid theory 

1. Cage size: order parameter 

The basic picture of the cloned liquid theorv |18l - l20j is to view the amorphous solid state such that each particle 
is frozen in a disordered configuration which looks exactly like a liquid but allowed to fiuctuate around the mean 
position within a finite range called as the cage size A. The cage size A is the order parameter of the theory which is 
determined self-consistently by extremizing a free-energy function G{A). The order parameter distinguishes between 
the liquid state A = 00 and the solid state ^ < 00. Such an idea has been proposed originally in the context of density 
functional theory for amorphous solids J85l - l88l | . 

An important feature of the order parameter A is that it can be viewed both dynamically and statically |87l . |88| 
much as the Edwards- Anderson order parameter of spin-glasses [8^ . First suppose that we monitor the dynamics of 
the system with the Hamiltonian Eq. Q keeping track of the positions of the particles Vi{t) [i ~ 1,2,..., N) as a 
function of time t. Then measuring the mean-squared displacement (MSD) of the particles, we may detect existence 
a metastable state if the MSD tends to saturate to a finite value in a large time limit, 

1 ^ 
A^ ^ hm — - 5]((r,(t) - r,;(0))2)i„i. (60) 

i—l 

Here (. . .)ini is the average over diffident initial conditions. Indeed the mode coupling theory (MCT) [Il|, II4I has 
detected a transition from the liquid phase ^d = 00 to a metastable solid phase ^d < 00 at the MCT critical 
temperature T^. Alternatively the same cage size can also be obtained in a static way using two replicas, say a and 
h, which have no mutual interactions but obey exactly the same Hamiltonian Eq. (j4]). 



^d=A^^|]((r«-r,^)2) (61) 

where r°(t) and Y\{t) are particle configuration of the replica a and b respectively. 

2. Mean-field scenario 

At the mean field level, the RFOT theory [a-S, Q assumes the following idealized features. First we assume that 
there are metastable states a = 1,2,... with free-energy N fa- The metastable states are characterized by the cage 
size A discussed above: 1) if we prepare the system in any one of the metastable states, the system remains there 
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such that Ad as defined in Eq. (|60)) is finite and 2) if we prepare two rephcas in that metastable state, they remain 
there such that Ag as defined in Eq. (|6ip is finite which takes the same value as A^. Then the total free-energy of the 
supercooled liquid state and the ideal glass state can be written formally as, 

- /3F(T) = log Z Z = Y, e-^^'f- = / d/e-^('^^-^(/'^» . (62) 

In the last equation the complexity (structural entropy) T,{f,t) is introduced as, 

iVS(/,i) = log^<5(/-/4r)). (63) 

a 

We assume that there are exponentially large number of metastable states in the supercooled and glass phase such 
that the complexity is finite. Apparently this description is the same as the free-energy landscape picture based on 
the notion of inherent structure [gj, ISa]- However, as we noted in sec. Ill A5[ the metastable state in the context of 
the mean-field theory corresponds to metabasins rather than individual inherent structures. 

At the dynamical transition T^, the complexity S > emerges. Usually the complexity !](/, T) is a function of / 
which is convex upward and vanishes continuously approaching some fnun{T) from above. (See Fig. [5]b)) The total 
free-energy of the system below Td is given by the saddle point value of the free-energy / which dominate the integral 
in the 2nd equation of Eq. (|62p . Similarly, the structural entropy is given by the value of the complexity at the saddle 
point. Upon lowering the temperature, the free-energy and the complexity at the saddle point decreases. Finally 
the structural entropy vanishes at the Kauzmann transition temperature Tk, i. .e the entropy crisis. To summarize, 
apparently different groups of metastable states become thermodynamically relevant at different temperatures in the 
supercooled liquid state T > Tk while the equilibrium statistical weight condenses onto a certain amorphous "ground 
state" at T < Tk. 

The scenario described above is known to hold precisely in a class of mean- field spin-glass models [5|, |^, ISO, IH^j Il24j 
in which the states a = 1, 2, . . . correspond to the solutions of the equation of states called as the Thouless- Anderson- 
Palmer equations. For the structural glass, the cloned liquid approach substantiated the scenario starting from 
microscopic Hamiltonian as we sketch below. 

3. Cloning 

Now let us recall the idea of cloning [13, l2J, |8J|- We consider m replicas (copies) a — 1, 2, . . . ,to of the original 
system which stay in the same metastable state (thus the name 'cloning'). The free-energy of the cloned system can 
be written formally as, 

-/3mF,„(r) = logZ„, Z,„ = ^e-"'^^/''(^)= f dfe-^'^^^f-^^f-^^l (64) 

Note that there is only one summation over the states a instead of m summations which is because we are assuming 
here that different replicas are staying in the same metastable state. We will recall shortly later (sec. IIIC Sp how to 
realize such a state in practice. 

Now we may regard ?n as a variational real number. Then we can notice that the saddle point fsp{T,m) which 
dominate the integral in the 2nd equation of Eq. (j64p can be moved around by varying m without changing the 
temperature T itself (See Fig. [5]b)). As we increase m, the saddle point /sp(r, to) becomes smaller and at some 
critical m*[T), the complexity at the saddle point Eq. (j65p vanishes, just like the physical {m = 1) Kauzmann 
transition. This observation suggests that even at very low temperatures below the putative Kauzmann transition 
temperature Tk, we may choose sufficiently small m such that the cloned system as a whole still remain effectively 
in a liquid phase where the structural entropy is still finite (See Fig. [5] a)). In the latter regime, we may compute the 
free-energy Fm (T) of the whole system regarding it simply as a liquid using the conventional liquid theory |23| as we 
recall shortly later. Then the value of the complexity can be evaluated as |2J], 

j:if,^iT,m),T) ^ - — ^— (65) 

The whole profile of the complexity curve J^{f,T) at a fixed temperature T can be "scanned" (See Fig. [5] b)) by 
varying m by which the position of the saddle point /sp(T, m), 

1 djmFrnjT)) ,_, 

UiT,m) = - g^ (66) 
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can be varied at will. Now we can locate the critical point m* (T) where the structural entropy vanishes 

I](/,p(r,m*),T)=0 (67) 

for each temperature T. Further increase of m does not change the saddle point value of the free-energy. This implies 
the value of the free-energy in the glass phase T < Tk is given by, 



lim F„,(T) = F„,.(T). 

The Kauzmann temperature Tk is located by (See Fig.[5]b)), 

m*(TK) = l 

At higher temperatures T > Tk, m* = 1. 

Thus the system in the glass phase T < Tk can be viewed as a liquid at an effective temperature, 



T* = 



T 



(68) 



(69) 



(70) 



m*{T) 

Typically it is found that ■m*{T) - T/Tk so that T* = T/m*{T) ^ Tk [H, El. This implies configuration of the 
system below Tk is almost like that of the liquid at Tk independently of the temperature T. On the other hand, in 
the supercooled liquid state T>Tk, m*{T) = 1 so that T* = T. 



a) 



b) 



S(/) slope m/T 




m* 1 



FIG. 5. Cloned liquid method: a) schematic m — T phase diagram b) complexity 



4- Disentangle intrastate and inter-state responses 



Let us make here a general remark that the cloning method allows disentanglement of generic linear response 
functions into intra-state and inter-state responses much as anticipated by the phenomenological free-energy landscape 
picture discussed in sec. Ill A 51 (See Eq. ((26)) ). In the present paper we fully makes use of this feature to analyze the 
rigidity of structural glasses. This is a strong advantage of the use of the replica method because the two responses 
are physically quite different in glassy systems: (3 and a relaxations as we discussed in sec. Ill A 41 In Appendix IBJ we 
demonstrate the method works exactly in a typical mean-field spin-glass model. 

Suppose that we have a physical observable of our interest, say O, which can be coupled to a conjugate external 
field h so that the Hamiltonian of the system can be written as, 



H = Ho- hO, 



l3FiT,h)=\ogZ Z = J2< 



-NPf^{T,h) 



(71) 



where Hf) is the original Hamiltonian of the system. Suppose that there are many metastable states a = 1, 2, . . . such 
that the free-energy of the system can be written as Eq. (jH)), which reads. 



(72) 



The static susceptibility. 



1 d{0) 



N dh 



1 d'^F{T,h) 



h=0 



N dh^ 



h=0 



- ^ iio') {O?) , 



(73) 
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can be formally decomposed as, 



with 



X^X + X 



x^IxJ x^^([o^]-[OoF) 



where we introduced averaging over different metastable states, 



...]^^P„(/.), 



Pc{h) 



'N0UiT,h) 



and defined thermal averages within metastable states, 

dUh) 



0« = {0)c 



N 



dh 



h=0 






(74) 
(75) 

(76) 
(77) 



ft,=0 



Here (. . .)« stands for a thermal average with respect to thermal fluctuations within the metastable state a. By 
inspecting the 1st and 2nd term in the decomposition of the susceptibility Eq. (|74p . namely x ^^d x which are 
defined in Eq. (j75|) . one can realize that they can be naturally interpreted as intra-state and inter-state susceptibilities 
respectively. 

Now let us turn to consider linear response of a cloned system. The Hamiltonian of the replicated system with 
replicas (a = 1, 2, ... , m) can be written as. 



B^ = Y, [{Hor - /i'O'^] 



(78) 



Here we have put probing field /i° on each replica (a = 1, 2, . . . , to). The free-energy of the cloned system reads, 

-/3TOF„(r,{/jn) =log^™(T,{/i"}) Z„(T, {/!''}) =^e-^^2::U/o(T,h") (^g^ 

a 

Note that different replicas stay in the same metastable state a but subjected to different probing fields h°' (a = 
1,2,..., to). Then we find a generic susceptibility of the cloned system, 



Xab 



md^F^{T,{h-}) 



N dh'^dhb 



{h-=a} 



^{{0^0')-{0^){0'))=XmSab + Xn 



with 



where 



IXa]T,„ 



Xrn ^ j^ {lOl]T,n - {O^l,., 



]r,m = ^-Pa(/i,r, to)... Paih,T,m) 



-Nml3fa,iT,h) 



Zm = J2 ' 



'NmPf^(T,h) 



(80) 



(81) 



(82) 



Here x a-nd x ^^e nothing but the intra-state and inter-state susceptibilities defined in Eq. (f75|) but at the effective 
temperature T jra. By taking trace over the replica index we find. 



X^Xab =Xm +TOXm- (83) 

b=l 

In the limit ra ^r \ this becomes nothing but the total susceptibility of the original system (see Eq. ([74])), 

=. lim Vxab = Xm-(r) +"^*(T)x™.(T)- (84) 

TTi— !-l 7 



X 



6=1 



Here we are assuming that we evaluate the terms on the r. h. s. in the liquid state (See Eq. 
Let us list below some observations and remarks. 
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• In the cloned "liquid" system, the total susceptibility Eq. (l84l) is just that of a (molecular) liquid state at 
temperature T* = T /m* . 

• By comparing Eq. (|80p and Eq. ([7^ we notice that replica trick allows us to disentangle the intra-state (/3- 
relaxation) and inter-state response (a-relaxation) hided in the (cloned) liquid: to find the intra-state suscepti- 
bility X one simply needs to look at the coefficient of the terms with 5ab in the susceptibility matrix Xab- 

• Note that the decomposition can be done not only in the glass phase where m = m*{T) < 1 but also at higher 
temperatures where m = m*{T) = 1 as long as the metastable states (with finite cage size A) exist. This is 
made possible by introducing m clones and taking to — > 1 limit afterwards (see Appendix |B]) . 



• In the standard interpretation of the RSB solution for the glass phase [90|, |91| , which is based on solvable mean- 
field spin-glass models ( see Appendix IBJ) . the factor to,*(T)(< 1)) which appears in front of x i^i Eq. (|M)) is 
understood as the probability that two independent replicas stay at different metastable states (here the replicas 
are not the 'cloned' ones) . In a more physical term, it roughly corresponds to the probability that a given system 
stay at an low lying 'excited' metastable state rather than 'ground' metastable state in the glass phase. The 
parameter m* (T) is decreasing function of the temperature T and typically behaves as ?7i ~ T JTy^ [92| . 

The idea presented above is generic but needs to be implemented more explicitly. In sec. Ill C 7[ we present a recipe 
to do realize the decomposition using the cage expansion. 

5. Molecular liquid 

Now we recall a microscopic implementation of the idea of the cloned liquid jlSf . The idea is to consider a molecular 
liquid state: a liquid with an assembly of hypothetical "molecules" i ~ 1,2, ... ,N with the center of mass (CM) at 
Ri, which itself consists of original particles of different replicas a = 1, 2, . . . , m at rf . The task goes as follows. First 
one put an artificial confining potential to a create molecular state: we start from a 'forced' cloned system whose 
free-energy Eq. (|M)) is given by, 

-. (. ^ ™ jd a 

- /?-^™ = l°g ^™ ^™ - n^iw / n n -^^-'"- (85) 

where the potential part of the Hamiltonian given by, 

in ^ N ni 

^»-EE-(4)-iEEir^-r?i'- (86) 

a=l i<j i=l a,b=l 

Here the 1st term on the r.h.s is the two-body interaction potential Eq. ^. The 2nd term is introduced to confine the 
particles in different replicas to develop a 'forced' molecular state. The cage size A of the resultant system is defined 
as, 

^^^^NdA ^Ij:± (« - r.^)^) ^ ^i^^^jfl (87) 

(=1 a, 6=1 ^ '^' 

After computing the free-energy Fm{a), a Legendre transformation, 

f3mGrn{A) = exta <^ l3mF„i(a) ^^ -Nd— 

[^ 2 a 

yields a free-energy function G{A). Here extr.Q stands for extrimization with respect to a. Subsequent extremization 
of the latter with respect to the cage size A, 



^ dGm{A) 



dA 



(89) 

^ — ^cxtr 



yields the cage size A = Aoxtr of the desired cloned liquid state, which remains after switching-off the artificial 
confining potential 1/a — >■ 0. Then the complexity can be computed using Eq. (|65p . The final step is to find m*{T) 
at which the complexity vanishes (See Eq. (p7)) V For clarity we recall the basic steps of the computation [l^, [l^ in 
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the following. Especially the parameter a — a{A) determined below as Eq. (|106p is indispensable in our analysis of 
the rigidity of glasses in sec. IIIII 

Anticipating that thermal fluctuations of the molecular liquid can be factorized into two parts: thermal fluctuations 
inside molecules and the thermal fluctuations of the CM positions of the molecules, we decompose the positions of 
the particles as, 

rf - R, + < (90) 

where u° are the fluctuations around the CM position of the molecules at R^. Then the integration variables can be 

changes as / Rti UZi d^rf = / Ul^ d^R. UZ^ d'^rfS''{R.-l, ^7=1 rf) = J Ul^ d^R, UZ, d^u^ruH^ (EZi <)■ 
The averages over thermal fluctuations inside molecules is given by, 

/N m , ,^ w„_]^s -d/2 / m \ 

n n '^'<^("i") ■ • ■ ^(>'^»^ (^^^^^) -5'^E""h"^^-^^^ '""""''' (91) 
i=la=l ^ ' \a=l J 

which yields in particular, 

{{u1Y{u]Y),^^, = -(1 - ni5ab)^5,j5^,. (92) 

Assuming first a is small so that the fluctuations m's are small one finds. 



-Nd 



1 / (2WA^)"'-^ ,N [^d'^R 
'" - TV! V 



/ n i\? n ^^p [-/^"i^cff (i?. -Rj)+..] (93) 



where we introduced the effective potential, 

t'eff(r) = v{r) - (1 - m)-^V\(r) + . . . (94) 

where V = {d/dx, d/dy, d/dz). This implies averages over thermal fluctuation of the CM positions of the molecules 
can be written as. 



V • ■)T/m,[v,!!(r)] 

with 



/ Yl dR,Wi{R,})T/m,[vMr)] ■ ■ ■ (95) 



e-/3i:.<,i'(fl.j) 



Now the free-energy of the molecular liquid can be expressed as. 



pmF,n{a) = TV log I y' ^^^"^^^T ' ^^ " PrnF,,^{T / m)[v,s{r)\ + . . . (97) 

Here the 1st term on the r.h.s is due to thermal fluctuations inside the molecules. The 2nd term Fna^iT / m)[v{ry\ is 
the free-energy associated with the configurations of the CM positions of the molecules which turns out to be just the 
free-energy of a liquid at temperature T* = T/m, 

-f3mFu^{T/m)[v{r)]^logZu^{T/m)[v{r)] Zn^{TMr)] ^ — /j] ^e"'^^-. ''(^-«^) (98) 

Assuming a is small we find. 



- /3m^nq(T/m)[veff(r)] = ~ iSmF^^iT / m)[v{r)] + /3(1 - m)-^ ^(V2«(i?,,)), + . . . (99) 

where we introduced a short hand notation, 

(■ • ■)* = (■ • ■)T/m,[v(r)] (100) 
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The resultant free-energy of the molecular state after the Legendrc transformation Eq. (|88|) is obtained up to 1st 
order in the expansion of the cage size A as, 



P ^ /,x «, / X dl — m) d(l— to), , ,,.9x ,1 — "^1 V^ ',/„9 / ^^ P ^ ,^ , , ^, ,9, 

^^^'"('^^ ^ " 2^;! ^°S(™) + 2m + 2to l°g(27rA/A2) - A— — - J^ /?(V2«(r,,))* + ;^^iiq(T/TO) + 0{A') 



i<j 

(101) 
More precisely the Legendrc transformation is done as follows: we expand the parameter a in power scries of A as 

[a, 

a = aiA + a2A'^ + a3A^ ... (102) 

The coefficients ai and 012 are fixed by the 1st order cage expansion while 013 requires the 2nd order cage expansion. 
Since we consider to, < 1, we have to maximize (rather than minimize) the free-energy GmiA) with respect to A. 
The equation = dGm{A)/dA Eq. ^ reads, 

= ^ - f ^ E^^'^'^'^'^))* + OiA). (103) 

Then we find the cage size of the metastable glassy state at the level of the 1st order cage expansion as. 

'^"2/3iE.<,(V2t;(r,,)),- ^^°^^ 

Finally the glass transition line to, = to*(T) (See Fig. [5|) is determined following the prescription discussed in sec. 
Ill C 31 To this end we use Fm{a{A)) obtained above in Eq. ((65)) and locate m = m*{T) where the complexity S 
vanishes. The Kauzmann transition temperature Tk is determined by Eq. (|69p. 

6. The effective Einstein model: thermal fluctuations of molecular liquid 

The analysis by the cage expansion implies the molecular liquid state can be described as follows. The coordinates 
of the particles are decomposed as Eq. (|90| which reads, 

where R^ is the CM position of the molecule and u,; is the deviation of the particle belonging to the a-th replica from 
the CM. The statistical weights of the configurations of the CM are given by Eq. ((M)) which is that of a liquid with 
renormalized potential Vcff{r) given by Eq. ([M)) . On the other hand, the statistical weights of the configuration of 
the displacements inside the molecules are given by Eq. (j9ip V Therefore the thermal average in the molecular liquid 
state becomes factorized as, 

/N m 

n d''R.WT,m,[vMr)] {{Ri\) n ^'<^({""}) (105) 

i=l a=l 

where (. . .)cagG is the average over thermal fluctuations inside the molecules defined in Eq. (|9ip ) and (. . .)T/m.[vcff(r)] 
is the average over thermal fluctuations of the CM positions of the molecules defined in Eq. (PSI) . 

The key parameter which defines the renormalized potential Vcsij") and the width of the fluctuations inside molecules 
is the parameter a. Going back to Eq. (|102p we find the parameter a = a{A) as, 



where the cage size A is given by Eq. (11041) . This follows from the fact that aiA = mA and a2A^ = mA (the latter 
can be seen using Eq.(38) and (40) of [ia|)- The result implies in particular that. 



{iufnu'^r)c.s. = {u')c.sA,d,^{m6ab~l) ("')cagc ^ -3^^^— — — -- = ^^^^ (107) 
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where we introduced 

Keff = —; (108) 

Eq. PITTI) follows from Eq. ^T^^Eq. ^ and Eq. piTO)) . 

It is instructive to note that the above result could have been obtained without using replicas but more intuitively. 
Suppose that we focus on a single particle, say i-th particle, and imagine that all other particles j{j^ i) are frozen in 
a certain configuration. Then the surrounding particles create a static effective potential Ui = ^ Si Mi) ^(''u) for the 
i-th particle: the Einstein model. Then the thermal average of the fluctuation of the particle can be evaluated as, 

knT 

Here we assumed isotropic fluctuations inside cages. Thus we find K^n defined in Eq. (|108|) can be viewed as the 
Hookian spring constant of the effective Einstein oscillator. 

The above discussion suggests that the basic idea of cage expansion is very similar in spirit to consider energy 
landscape in the vicinity of inherent structures (see sec. Ill ASp . However there is an important difference : the 
configurations of the center of mass positions of the molecules Rt are by no means restricted to energy minima but 
allowed to take any configurations of a liquid thermalized at effective temperature Toff = T/m* (see Eq. ([95]) ). As we 
discuss soon, this point is important in linear responses to external perturbations at finite temperatures. 

Using the statistical weights Eq. (|105p thermal averages of physical quantities can be done as follow. Consider 
an observable 0{{rfA) which is a function of the distances between particles r°- = |r" — r°|. Using the molecular 
coordinate Eq. (|90p and assuming that the fluctuations uf's are small, it can be expanded as, 
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i<],k<l a,b ^ ^i^ijY^i'^kl) 

Then the thermal average can be done as. 



\ E T:£ ^!^,M -^r^-^'-- im 



l'2x Y- Y- V- Y-n . ^Y-/ ^'^(K,}) 



{0{{rf,}))^ (0({%}))t/™j..,,M] - ^(«')cage }^}^ E E(l - ^^^-b) 



= {0{{R.,}))* + /3*(1 - m){u'U,,Yl m{R^J})yMRu))* {0{{R^,}))A^MR.,))A +■■■ 



' Jl(#j)j2(#i) a,b 



2^-.^a,.^^ ^.^- —'^\d(rf^^Yd(rt^J, 



where a/m? = 2A/m as given in Eg. (|106p . The evaluation of the averages (. . .)* (see Eq. (|100p ) can be done using 
methods of standard liquid theory [23| at temperature T/m. The value of the cage size A = A{T) and the parameter 
m ~ m*{T) must be determined beforehand. 

7. Evaluations of intrastate and interstate thermal fluctuations by the cage expansion 

In sec. Ill C 41 we discussed the possibility to disentangle the intra-state and inter-state responses using a cloned 
system. Here we present a more precise prescription to do it explicitly using the cage expansion. In Appendix |B] we 
demonstrate that the prescription works exactly in a mean-field spin-glass model. In the next section we compute the 
rigidity applying the following prescription to evaluate shear-stress correlation function which appears in the static 
fluctuation formula of the rigidity Eq. ([8]). 

Here let us consider the generic susceptibility considered in sec. Ill C 41 (see Eq. ([501)). 

Xab=J^{{OaOb)~{Oa){Ob)). (112) 
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Using the above prescription Eq. (|llip it becomes, 



Xab = -^ ((0^({i?,j}))T/mJu„ff(r)] ^ {Oi{Rij}))T/rn,[v,t{ir)]) 



' N 



EE 



do do 

d{xlY d{x']Y 



{l-m6ah){{u;r{ujY) 



31 /cage 



(113) 



with {{ui)'^{ujY) cage given by Eq. (|107p . Comparing the above resuh with the last equation of Eq. (|80p we find the 
intra-state susceptibihty as, 



X- 



P* 



EE 



do do 



d{xlY d{x]Y 



{{u^nu,r) 



(5* 



EE 



do 



with (3* = m/3, and the inter-state susceptibihty. 



(114) 



X = I ({0'i{R^,}))T/m,Mr)] ' {Oi{R^3}))T/rn,Mr)]) " |EE 



do 



I cage • 



(115) 



FinaUy the total susceptibility Eq. ([84)) becomes 

/3* 



X = ^ ((0^({i?^j}))T/m,Kff(r)] - (0({i?ij}))T/,„,[«„„(0] j 



(116) 



Although we have only performed the cage expansion up to 1st order in Eq. (|114p .Eq. (|115p and Eq. (|116p . we 
expect that the basic structure of the result is generic. Higher order cage expansion will yield more accurate effective 
potential v^sir) [201 which will refine the averages over the CM positions of the molecules. 



8. Discussion: responses inside metastable states 

Let us focus on the responses within metastable states and discuss it more in detail. Since fluctuations in the close 
vicinity of an IS is completely characterized by the elastic hamiltonian Eq. p4p . the susceptibility xig associated with 
an IS can be obtained just like the 1st equation of Eq. pi4p but replacing the correlator {{uiY {ujY) cage by, 



{{u.T{^ 



ns 



{(3H- 



)ij 



with H ^ being the inverse of the Hessian matrix Eq. ([55]) . As the result we find, 



»S = T7EE 



do 



N '^^-' dix'^Y 

ij fj.u ^ '' ' 



do 



IS 



9(x5)- 



(ii-X 



(117) 



(118) 



IS 



where the derivatives are took at the IS. This must be compared with the intra-state susceptibility x defined in 
Eq. (|114p together with Eq. \\(Y1\ . which reads, 



N 



EE 



do 

dixfY 



1 

Koff 



(119) 



where we used Eq. (|108p . 

Comparing the two susceptibilities xis s-nd x ^^ find that they are apparently similar to some extent. In both 
cases, the fiuctuation of the observable O is induced by small fluctuation of particles around some reference state(s). 
We also notice that the cage expansion approach amounts to replace the full Hessian matrix by that of an effective 
Einstein model we discussed in sec. Ill C 61 Validity of such an approximation should be examined critically. A better 
approximation may be obtained for example by doing partial harmonic re-summation exactly |18| . 

However here let us focus on a more important difference between the two susceptibilities that comes from the 
difference of their reference state(s) themselves. Note that the reference state for x ^-s defined in Eq. pi9p is a 
thermalized liquid configuration within a metastable state at a finite temperature T* = T/m*{T), where all sorts of 
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thermal excitations, not only the elastic but also plastic deformations, are included. On the other hand the reference 
state for xis is just an energy minimum. This implies for instance, 



do 



d{xty 



IS 



< 



' do 



(120) 



The above observation suggests that in general the intra-state susceptibility x is larger than that of the inherent 
structures xiSi 



X > Xis- 



(121) 



In sec III A 5] we argued that a mctastable would be identified as a metabasin (MB), which is a union of inherent 
structures (IS), rather than a single IS. Then the difference between the two susceptibilities Eq. (J12ip should be 
attributed to thermal fluctuations between different inherent structures within a common metabasin, which are 
necessarily plastic deformations. In the supercooled liquid state T > Tk, m*(T) = 1 so that T* = T. This implies 
the difference between the two susceptibilities increases with temperature at T > Tk- On the other hand, in the glass 
phase T <Tk one typically finds m*{T) - T/Tk so that T* = T/m*{T) - Tk 0,[il. This suggests that x becomes 
nearly independent of temperatures below Tk such that x ~ Xis- 

We argued in sec. Ill A 41 that the same intra-state susceptibility can be observed either in equilibrium or out-of 
equilibrium based on the observation that in general the /3-relaxation shows little aging effects. In this respect, it is 
interesting to note that the intra-state susceptibility x reflect the underlying equilibrium state. Especially the above 
discussion suggests even the Kauzmann temperature Tk itself may be reflected in the temperature dependence of %. 



III. RIGIDITY OF GLASSES 



A. Shear on the cloned liquid 



In sec. Ill A H and sec. Ill A 21 we discussed some general aspects of static response to shear. The rigidity is defined 
as the 2nd derivative of the free-energy with respect to the shear-strain 7 as in Eq. ([2]) whose explicit expression, the 
fluctuation formula, is given in Eq. (j8]). Now let us consider to shear the cloned liquid. Suppose that we have a cloned 
liquid of m replicas which are subjected different shear-strains 7° (a = 1, 2, ... , m): the replicas are put into different 
containers which can be deformed independently. This is just an example of the static perturbations on the cloned 
liquid discussed in sec. Ill C4I Then we may formally expand the free-energy of the cloned liquid as. 



N 



mF™(7) = mFraiO) + iV^(a'')7, + y ^ MabT^T' + 0(7'') 



(122) 



;,6=1 



The argument for the decomposition of the response Eq. (|80p implies that the generalized rigidity fiat may also be 
parametrized as. 



^^ah = iiSa,b + P- 



with the intra-state rigidity or effective rigidity, 



A* 



5^/0(^,7) 



c}-y2 



7=0 



and the negative correction term due to inter-state fluctuations, 
where a^ is the 'internal stress' of state a, 

a/a(T,7) 



97 



(123) 



(124) 



(125) 



(126) 



7=0 



As we discuss below fi, which is negative, can be interpreted as inter-state modulus. 
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F(7) 




FIG. 6. Schematic mean-field picture of free-energy landscape and static stress-strain curve below Tk. The metastable states 
(metabasins) with rigidity jl are represented by parabolas 



The trace over the replica index of pLab yields the rigidity of a molecular liquid. This must be zero since cloned 
liquid as a whole is just a liquid (!) (see the discussion in the end of sec. Ill A2[) 



^ Mab = A + "^M = 0- (127) 

6=1 

Following the discussion in sec. Ill C 41 we find the total rigidity of the glass from Eq. (|127p as, 

/i = lim y ^ab = (1- + ""Ti* fi = Q- (128) 



b=l 



where m* = m*{T) is the glass transition line. Thus we find the inter-state correction term fi is simply related to the 
intra-state rigidity (i as, 



1 
m* 



(129) 



B. Static analogue of yielding or shocks 



Let us discuss a physical interpretation of the above result Ecj^. p29p . Suppose that internal stress a a defined in 
Eq. (|126p has a Gaussian distribution with zero mean and variance A'^ /N. Then Eq. (|125l) and Eq. (|129l) implies that 
the intra-state rigidity is related to A asJ49f. 



fi^m*/3A^. 



(130) 



Note that this corresponds to Eq. ([5^ . (We mentioned similar relations for t he in stantaneous rigidity Eq. P7)) and 
the rigidity of inherent structures Eq. pip ). A very recent numerical study 112| explicitly demonstrated that the 
internal stresses of inherent structures indeed follow Gaussian distributions. 

The static response of a system with many metastable states against variation of the shear-strain 7 at low temper- 
atures T < Tk niay be viewed pictorially as shown in Fig. [51 Variation of the free-energy fa{"f) associated with a 
given metastable state a along the 7-axis may be described as. 



/a(7) = f(7-7r')' + r 



(131) 



where /t is the intra-state rigidity. The free-energy takes a minimum value /™'" at 7 = 7™'"- Then the internal stress 
fa (7) of the metastable state at a given 7 is given by, 



ao(7)=A(7-7r") 



(132) 



Of course wc should expect that the metastable state itself disappears at some distance from 7™™ approaching its 
metastability limit where a saddle- node bifurcation takes place [63 1. 
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Let us consider a set of metastable states which lives at around a given point 7 along the 7-axis. The locations 
of their minima 7™™ will be distributed randomly around the point 7. Then the internal stress (Jail) defined by 
Eq. (|132|) at a given 7 will take both positive and negative random values depending on the sign of 7 — 7™'"- Then 
the relation Eq. (|130|) implies the variance of the distribution of the distances 7 — 7™™ is {\/m* l3jl)N~''^ . 

At any point along the 7-axis, there will be exponentially large number of metastable states oc e^^^^' at a given 
free-energy level (per particle) / where S(/) is the complexity or structural entropy Eq. (j63p . However, at low 
temperatures T < Tk, the ground state and at most only a few other low lying states would have significant statistical 
weights |92| . Under variation of 7, there will be level crossings among such low lying states as shown in Fig. [51 
Just around the points where there are level crossings, the rigidity should become locally negative. Correspondingly 
the ground state energy exhibits a cusp singularity where a level crossing takes place. This would be the physical 
interpretation of the negative inter-state modulus p, ~ — (l/m*)/i (Sec Eq. (|129p ). We can consider this feature as a 
static analogue of yielding. 

Typical spacing between such level crossings along the 7-axis will be of order, 

'-'■ - ^- ('-' 

This is because typical free-energy gap between low lying metastable states will be of order /cbJk while free-energy 
change of each of them for a small change of the strain 7 is proportional to their internal stress CTq, . Another important 
scale is typical width of the thermal rounding of the level crossings. This would be simply give by. 

In the zero temperature limit T — > 0, the free-energy landscape should exhibit cusp singularity. Correspondingly we 
find that the intra-state modulus jl defined in Eq. (I129|) diverges negatively because limT^o rn*{T) — >■ 0. 

Based on the above observations the probability ^crossing to meet a level crossing at a given 7 may be estimated as 

^crossing ^ 5-f^ / hs ^ T /Tk (135) 

It is interesting to recall that the factor m* (T) which appears in front of the inter-state modulus jl in Eq. ()128p can 
be regarded as the probability to make an inter-state transitions within the standard interpretation of the replica 
symmetry breaking ansatz [90,|91|. (See sec. Ill C 41 and Appendix [B] for related discussions) Typically it is found that 
m*{T) ^ T/Tk [la, 1131- Thus Eq. (|135|) is consistent with the standard interpretation of the RSB ansatz. 

Note that in the thermodynamic limit N -^ oo, the spacing between the cusps vanishes [9J, [OJl- It means that if 
we define a coarse-grained rigidity for any small but non-zero interval along the 7-axis, we obtain rigidity in the 
iV — >■ (X) limit. Thus the static intermittcncy of the stress-strain curve sketched in Fig. [6] explains how we can solve the 
apparent paradox discussed in sec. Ill A II in the static mean-field theoretical context. Similar static intermittencies 
have been found in disordered systems including spin- glasses J93l - l98j . and elastic manifolds in random media 99l - ll01 |. 
which are actually related to the problem of intermittcncy in turbulent flows due to shocks [102| . 

The above observation implies the intcr-state fluctuations can be considered as Goldstone modes characteristic in 
glasses which try to restore the broken translational symmetry. On the other hand, the elastic fluctuation within a 
metastable state is more usual Goldstone mode which is also present in crystals. 

Although we have assumed T < Tk in the above discussions, the transitions between different metastable states 
under external strain can be relevant even at higher temperatures in the following sense. Within the so called mosaic 
picture of the RFOT theory [3, |^, lla, [131 , the conflguration of the system in the intermediate temperature range 
Tk < T < Tc consists of mosaics or patches of local regions such that the subsystem in each of the patch stays in a 
metastable state. The size of the local regions is called as the mosaic length ^mosaic- It means that the subsystem 
within a mosaic is as if at a temperature below Tk. 

Then an interesting possibility is that the intermittent stress-strain curves shown schematically in Fig. |6] may be 
observable by some rheological experiment in a confined geometry. Suppose to exert some local shear-strain 7 in 
a quasi-static manner on a system confined in a narrow region of size L. If the size of the region is small enough 
such that L < ^mosaic, then the resultant stress-strain curve may exhibit the intermittent feature with typical spacing 
between the shocks ^7s oc 1/vZ^. 

Finally also note for clarity that yielding at mesoscopic scales can occur by the saddle-node bifurcation mechanism 
at metastability limits [32[, [3^, [Z3|- I^ ^^^ energy level at the metastability limit is order 0(1), we find again that 
(57s (X l//z7. 
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C. Cage expansion of the rigidity 

Now let us examine the rigidity fiat defined in Eq. (|122p . Its microscopic expression can be obtained generalizing 
the fomulac Eqs. Eq. (IHl),Eq. © and Eq. ^ . 

d'^mFra 

^^<^>> ^ d-y'^djb "" %'--}])'5a6 - ^[(cr[{r»}]Cr[{r^}]) " (cr[{r-}] ) (cr[{rj}])] • (136) 

where the subscript, for instance, [{rf}] indicates that coordinates r^ (i == 1,2, . . . ,N) of particles in replica a. The 
1st term b on the r.h.s of the above equation is the born term defined in Eq. ^ which reads, 

h{r.}] == ^E^(^^^) ^(^^^■) = 4 [r'v^'Hr)S:l+rzM\r)il~£l)]^^^ . (137) 

The 2nd term is the correlation function of the fluctuation of the shear-stress, 

'^[{r^}] = irf^'^i^^j) cr(ry) = ru(^)(r) i,j%. (138) 

N ^ 

Our task is to evaluate the rigidity fiab by the cage expansion as follows. It is essentially the same as the analysis of 
the generic susceptibility Xab we discussed in sec. Ill U 71 First we explicitly decompose the coordinates of the particles 
as Eq. daOl), 

Then we expand the thermal averages which appear in the fluctuation formula of the rigidity Eq. p36p in power series 
of the molecular coordinates uf (a = 1, 2, . . . , m, i = 1, 2, . . . , N) using the formula Eq. (jllip . Then we plug-in the 
expansions in the microscopic expression of the the rigidity fiat given by Eq. ()136p . Finally we evaluate the thermal 
averages using Eq. ([9T|) (Eq. (|92|) ) together with Eq. (|106p for the molecular coordinates u° and Eq. (jlOOp for the 
center of mass coordinates R^ of the molecules. The result may be put into a power series of the parameter a (or the 
cage size A), 

llab^ {lJ'0)ab + {iJ^l)aba + O{a^) (139) 

or 

fJ-ab ^ fiSab + fj- fi = fiQ + flia + 0{a'^) (1 ^ jlo + jlia + 0{a'^) (140) 

which must satisfy the sum rule Eq. (|127l) . 

^{^J■n)ab = An + rnp.„ = (141) 

b 

at each other in the cage expansion. 

1. Zero-th order 
At zero-th order of the cage expansion, all particles are bound to the CM positions of the molecules. We find, 

itJ-o)ab = {b[{Il,}])Jab - I3[{<J[{R,}])* - (ct[{r,}]>*] (142) 

which obviously satisfies the sum rule Eq. (|14ip . 

m 

X!(/^o)a6 = (%R,}])* - I3m[{af^^^y]),. - (ct[{rj])J] = 0, (143) 

a=l 

since this is just the rigidity of a simple liquid at T = T/m*. Thus at the zero-th order we find the intra-statc and 
inter-state rigidities as, 

Ao == (^[{R,}])* Ao = -■ (144) 
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2.. First order 

At the 1st order of the cage expansion we find the following result, 

{^J■l)aba = Ji(l - mSab) + (1 - m)[{J2 + ,h)5ab + Ji + Jb] 



(145) 



where 



•^1 = ^^ E E E /3 (V<T(r,,, ) • Va(r.,, )) , 



* ji(7^*)i2(#») 



^2 = -^^E(^'^(r^^))* 

i<j k<l 

n<ji i2<J2 
■^5 = -'^^ E E E/3/5*[(a(r,,,Ja(r.,,JV2^;(r,0)* - (a(r,,, Ja(r,,, J).(V2z;(rfeO>*] (146) 



N 

n<ji «2<i2 *:<; 

Here we introduced a parameter, 



C= (w^)c 



;3*K, 



cff 



(147) 



with Keff defined in Eq. (|108p . In the above equations Eq. (|146p and in the following quantities which appear inside 
the brackets (...)* are functions of only the CM coordinates R,; (i = 1, 2, . . . , N). 

As discussed in Appendix [Dl the group of terms J2,J3,J4,J5 can be related to the rigidity of renormalized liquid 
interacting with the effective potential Uoff(?') given in Eq. (|94p . Then the following identity. 



J2 + J3 + m(J4 + Js) =0 
holds. Using this identity, the result Eq. (|145p can be simplified as, 

1 — 771 



yLia = {jliSab + fj-)a 
so that 

fiia ~ — m 
The result satisfies the sum rule Eq. p4ip 



Ji - 



m 



{J2 + J3) 



(1 - mdab)- 






m 



^{fJ'l)ab = 



as it should. To summarize, up to 1st order in the cage expansion we find. 



(148) 



(149) 



(150) 



(151) 



(152) 



3. Summary and discussions 



To summarize, we find the rigidity of the cloned liquid Eq. p36p as. 



fJ-ab = K^ab ) 

777 



(153) 
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FIG. 7. Fluctuation of stresses due to fluctuations of particles inside a cage 

In practice, it is useful to express various terms in Eg. (|152p in terms of particle distributions so that standard 
techniques of the density functional theories of liquids j23| can be incorporated. The details of the technical aspects 
are presented in Appendix |e1 [F] and [G] 

Let us discuss the two terms on the r. h. s. of Eq. (|153|) below, 

• Born term of the rcnormalized liquid 

The terms J2 and J3 are corrections to the Born term due to the renormalization of the potential Eq. ([94]). It 
is easy to realize that the sum of the three terms, i. e. the Born term, J2 and J3 is equivalent to the Born term 
of the renormalizcd liquid with the effective potential Eq. (|M)) . 

(&[{R.}])T/m,K«(r)] = (fe[{R.}])* + (1 - m){J2 + J3) + 0{A^) (154) 

with (. . ^T/mAv^atr)] defined in Eq. (j95p . Here ^ff^^^ii is the Born term of the renormalized liquid which can be 
obtained by substituting Woff('') into v{r) of Eq. ^. 

• Non-affine correction 

The term Ji is clearly distinct from others. While the contributions of the terms J2 and J3 vanish above Tk 
where m*{T) = 1, the contribution of the term Ji docs not. Presumably it continues to exist up to T^ where the 
metastable solid states disappear. Note that it has the same structure as the generic intra-state susceptibility 
X (See Eq. (|114p ) discussed in sec. Ill C 71 Let us rewrite the Ji term as. 



^^ = 7^^E(l-'a (155) 



Kcff N 

where KcS is defined in Eq. (jlOSp and H is the derivative of the shear-stress introduced in Eq. (|38p . In [63] it is 
pointed out that the field H^ defined in Eq. psp is a random field with very short-ranged correlation in space. 

The physical interpretation of this term is straightforward and quite suggestive. Consider time scales shorter 
than the a- relaxation time, where the stress- field does not evolve much so that, as a first approximation, we can 
assume it is a frozen-in vector field associated with a metastable state. As shown schematically in Fig. [71 the 
fluctuation of a particle, say particle i, inside its cage induces fluctuations of the stresses between the particle 
and the surrounding particles, say j'l and j2- The point is that the resultant fluctuations of the stresses aij-^ and 
Gij^ are correlated since they originate from the common source of fluctuation, the fluctuation of particle i. 

Although we have only performed the cage expansion up to 1st order, we expect that the structure Eq. (jl53p is 
generic. Higher order cage expansion will yield more accurate effective potential Vcji{r) |20|. This will amount to 
reflne the renormalizcd Born term Eq. (|154p and the Ji term. Similarly, the S fleld Eq. psp will be replaced by that 
of renormalizcd potential at higher orders in the cage expansion. 

It is instructive to compare the above result with that of the rigidity /iis of inherent structures (IS) discussed in 
sec. Ill A 71 The comparison is essentially the same as the comparison of the susceptibility associated with inherent 
structures xis and the intra-state susceptibility x ^^ discussed iu lIICSI Comparing the full expression of the rigidity 
of inherent structures /xis Eq. ([37]) with our results of the intra-state rigidity /t Eq. (|155p together with Eq. (|147p . 
we flnd that the our scheme approximate the full Hessian matrix by that of an effective Einstein model. The other, 
more important difference is that the reference state itself used in the evaluation of the intra-state susceptibility fi is 
not completely frozen-in one of energy minima, at marked variance with that for ^i^. The reference state is just took 
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from a thcrmalized liquid configuration in a mctastablc state. Then it is natural to expect that the amplitude of the S 
field (I Hip) , which appears in Eq. (|155p . is larger than that at inherent structures (See Eq. (|40p ). Thus we consider 
that the intra-state rigidity fl obtained above is most likely smaller than the rigidity of inherent structure /ijs for the 
same reason as the generic intra-state susceptibility x is larger than the susceptibility xis associated with inherent 
structures as we discussed in sec. Ill C 81 We consider that the underlying physical processes are thermal excitations 
between different inherent structures within a common metabasin, which are probably localized plastic deformations 
as we discussed in sec. Ill A 81 

Finally let us recall the discussion in sec IIIBI where we pointed out that mean-field theories generically imply that 
the rigidity of solids is closely related to the order parameter, more precisely the rigidity is simply proportional to the 
square of the order parameter (see also Appendix. [C] for a corresponding problem in a certain vectorial spin-glasses) 
. Indeed the computation of the intra-state rigidity fl of glasses performed in the present section is consistent to some 
extent with this generic picture. We found the intra-state rigidity fl can be expanded in power series of the order 
parameter A as Eq. (|140p combined with Eq. (|102p ). which reads, 

fi{A)==co-ciA + 0{A'^). (156) 

with Co and Ci being positive constants related to /tg and —fti respectively (See Eq. (|144p and Eq. (jlSOp ). The 
above expression Eq. (jl56p implies larger cage size A yields smaller rigidity fl. Here it is instructive to note that the 
equivalent of the cage size A in the spin-glass case is 1 — g with q being the Edwards- Anderson order parameter as 
we point out in Appendix IB] (see also Appendix IC)) . 

D. Fluctuation of the shape of the container : shear-strain fluctuations 

So far we have implicitly assumed that the shape of the container, parametrized by the shear-strain 7, is rigid and 
fixed. What will happen if we allow it to fluctuate? As long as the rigidity fi is positive, we can safely integrate over 
the shear-strain 7 to obtain an elastic contribution to the free-energy of the cloned liquid. Using Eq. (|153p we find. 



N J a V «'""' (/3A(a)/H'"" 



/Xab = fiia) [Sal, ■ (158) 



with 



Here the intra-state rigidity fi is given by Eq. (|140p . 

fi{a) = fio + fiia + 0{a'^). (159) 

1. First order 

Let US first discuss the consequence of the fluctuation of the shear-strain 7 at the 1st order in the cage expansion. 
By adding F^^^^^'^{a) to the free-energy Fm{a) Eq. (jHZ]) and performing the Legendre transform Eq. ([55)1 we find 

a{A) = mA+ | ^1 ^(vM^..))* + ^| J A' + 0{A^) (160) 



Then the free-energy Gm{A) Eq. (jlOip is modified as, 

4g,„(A) = t^ log(2.A/A2) - aI^^I-Y: P{yMn,))* V^ logimc^iA))) + . . . , (161) 

N 2m m N ^-^, 2m 



where we omitted terms independent of the cage size A. Now the equation = dGm{-A)ldA by which the cage size 
is determined Eq. (|103p is replaced by, 

°4-f^S<VM..,».-|f + OM (162) 

l<3 
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which yields, 

^ ^ (163) 

2/3wE.<,(V2«(r,,)>*+™Ai/Ao' 

Since fii is negative (see Eq. (|150p ) the result Eq. (|163p means that the fluctuation of the container tends to increase 
the cage size and thus weakens the glassy order. 

In turn, the increase of the cage size will make the rigidity smaller. Using Eq. (|163p back in Eq. (|160p . we find 
again a = 2mA, which is the same as Eq. (|106p . Thus the intra-state rigidity of the system which is contained in the 
deformable container becomes, 

fi — fiQ — 2mAfli. (164) 

This is formally the same as the previous one Eq. (|153p but it is now smaller than the latter because the cage size A 
is increased as found above Eq. (|163p . 

2. Coupling of the cage size and the shear-strain fluctuations : a possible scenario 

An interesting feature found above is that fluctuation of the shear-strain is coupled to the cage size A, which is 
the order parameter of the glassy state, and makes it larger. We suspect that the actual effect of the shear-strain 
fluctuation would be much stronger than the behaviour found at the 1st order level. Let us discuss a possible 
consequence of the coupling below. 

By fully taking into account the elastic free-energy beyond the 1st order cage expansion, the saddle point equation 
for the cage size A Eq. (|162p may be cast into the following simple form, 

11 r 1 

= T--r + ^-r^ (165) 

A Aa l-r^Ao' 

where Aq is the cage size in the absence of the shear-strain fluctuation. We also introduced a dimensionless parameter 
r which represents the strength of the non-afhne correction, 

r=-^m^o- (166) 

Mo 

Note that /io > and /ii < (see Eq. (|144p and Eq. (|150p ) so that r is a positive parameter. We anticipate that it is 
an increasing function of the temperature T. Although we have not performed higher order cage expansions explicitly, 
we expect that the structure of Eq. (|165p is generic. It can be simplified as. 



where 



The solution of Eq. p67p is easily found as. 



(167) 



a^A. (168) 



A 







vT^4f , , _ ^, 2 



2r 



l + r + 0{r^) (169) 



where we have took the physical solution which becomes 1 in the limit r — > 0. Most interesting feature is that the 
cage size exhibits a square- root singularity approaching a critical point. 



a-«c(l-J^-l) (r^r-), (170) 



with Tc = 1/4 and Oc = 2 and disappears at Tc- Correspondingly this is reflected on the intra-state rigidity as. 



/i~/ic(l + jf -1) [r^r-), (171) 
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where jlc = /io/2- Thus the rigidity exhibits a discontinuous jump preceded by a square-root singularity. 

The above observations imply the consequence of the coupling between the order parameter and the strain fluctua- 
tion can be significant leading to collapse of metastable solid states when the strength of non-affinity r become strong 
enough. Interestingly enough the the discontinuous vanishment of the rigidity preceded by the square-root behaviour 
Eq. (|f 7f I) is apparently similar to the expectation in the pure mean-field limit discussed in sec. IIIBI However here 
the rigidity is not merely subjected to the order parameter at variance to the pure mean-field scenario but plays an 
essential role for the melting, partly supporting the intuition of Bern's rigidity crisis scenario (37| . 

E. Discussion: how the rigidity disappears? 

An interesting general question is how the rigidity of the metastable solids disappear by raising the temperature. Of 
course glasses melt at the glass transition temperature Tg. But the effective rigidity fi remains finite in the supercooled 
liquid state above Tg so that a supercooled liquid can be viewed as a solid with finite life time. Our interest is how 
the rigidity of this temporally solid disappear by increasing temperature. 

1. Mean-field theories 

In sec. IIIBI we pointed out that mean-field theories generically imply that the rigidify of solid states are simply 
proportional to the square of the order parameters. Then it immediately follows that the rigidity vanish discontinuously 
signaled by preceding square-root behaviours approaching spinodal like melting temperatures from below, simply 
reflecting the same behaviour of the order parameters. We showed that this scenario holds in a Ginzburg-Landau 
theory for superheated ferromagnet (sec. IIIB"T|) . the mode coupling theory for supercooled liquids (sec. IIIB"^ and a 
replica field theory for vectorial spin-glasses (Appendix. [C]) . 

In sec unci we performed a microscopic computation of the rigidity of supercooled liquids and glasses based on the 
cloned liquid theory. Indeed the result also suggests that the rigidity is intimately related to the order parameter A 
as Eq. (|I56|) . which reads as, 

fi{A)^cn-ciA + 0{A^). (172) 

with Co and cjbeing certain positive constants. Since both the mode coupling theory |lll . |l2| and the cloned 
liquid theory [20[ predict that the cage size A exhibits a discontinuous jump preceded by a square-root behaviour 
A{Tc) — A{T) ex \/Tc — T as T ^ T~ , the rigidity fi as given by Eq. (|172p should also exhibit a discontinuous jump 
at Tc preceded by a square root singularity, at least at the mean-fie ld le vel. Let us note that a recent formulation of 
a replica approach to the rigidity of glasses by Szamel and Flenner |103l | is consistent with this scenario. 

To summarize, at the pure mean-field level, the rigidity does behave as the order parameter. Said differently, at this 
level the transverse fluctuations such as shear-strains, spin-waves, e.t.c. do not play any significant roles compared 
with the longitudinal fluctuations such as the fluctuation of the density, amplitude of vectorial order parameters, e.t.c. 
because the rigidity is totally subjected to the order parameter which parametrize the strength of the longitudinal 
fluctuations. The situation may change if we consider explicitly influences of the transverse fluctuations going beyond 
the pure mean- field level as we discuss below. 

2. Influence of transverse fluctuations 

In sec IIIIDl we discussed an alternative melting mechanism due to coupling between the order parameter, i. e. 
the cage size A and fluctuations of the shape of the container, i. e. shear-strain. Interestingly enough it also 
predicts discontinuous jumps of the cage size A and the rigidity fi with square-root singularities. However the driving 
mechanism of the melting is different from the one discussed above. In this case the control parameter is the non- 
affinity r defined in Eq. (jI66p which is presumably an increasing function of the temperature T. 

At the pure mean-field level, fluctuation of the shear-strain can be forgotten by just considering a container with a 
fixed shape. However we may view a finite dimensional system as an assembly of deformable cells for each of which 
we can define an effective mean-field free-energy. Such a local mean-field free energy on a cell around a point r may 
be expressed in terms of the local order parameter A{r) and the shape of the cell parametrized by local shear-strain 
7(r). In this description the shear-strain 7(r) naturally becomes a fluctuating continuous fleld (elastic strain fleld) 
which wc must integrate out. 

The equivalent of the shear-strain field 7(r) in spin systems are the derivative of the angular field V0(r) which 
parametrizes the spin- waves discussed in sec III B T] (fcrromagnets) and Appendix [Cj f spin-glasses). Apparently the 
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the spin-wave V0(r) must be integrated out in order to take into account fluctuations around the mean-field theories. 
As the result we will obtain certain renormalized effective potential for the longitudinal fluctuations, which should be 
analogous to the one we discussed in sec. IIII Dl 

However, as long as the effects of the transverse fluctuations can be took into account perturbatively, the overall 
basic characters of the spinodal like criticality themselves would remain intact. It will be very interesting to examine 
possible non-perturbative effects of the transverse fluctuations in future. 

3. Comparison with related problems 

The spinodal like criticality is very different from the case of usual 2nd order phase transitions (e.g. 0^ theory) 
where the order parameter behaves continuously at the transition point so that the rigidity also exhibits a continuous 
behaviour across the transition temperature as we noted in the end of sec. IIIB II In the latter case the vanishing of 
the rigidity is related to the divergence of static correlation length at the transition through the so called Josephson's 
scaling relation [52| . Such co ntinu ous behavior of the order parameter and the rigidity have been predicted for the 
cases of spin-glass transitions jl04| and vulcanization transitions |44| - t46| . 

Jamming transition is quite intriguing in this respect because it cannot be understood simply along the line of 
thoughts discussed above: some simplified model systems for granular matters exhibit discontinuous jump of the 
contact numb er preced ed by a square-root singularity but continuous behaviour of the rigidity at the (un) jamming 
transition H, E^, [lOi . 



F. A model computation - binary soft sphere system 

We now perfo rm an explicit computation of the rigidity for the case of a binary mixture of soft-spheres in three 
dimensions |106| . which is a standard model system to study the supercooled liquid state and the glassy state. We 
compare our results with the rigidity of the same system obtained by Barrat et. al. [36| using molecular dynamic 
(MD) simulations. Thermodynamic properties of the same system has been analyzed using the cloned liquid theory 
by Coluzzi et. al. ^. 

1. Model 

The system consists of two kinds of particles e = -|-, — which have different radius CT-j. and ct_ and interact with 
each other by a repulsive potential, 

.^^^.(.,,.)=«„(^^!i±^y'. (173) 

We introduce a smooth long-range cut-off of the potential following |43| which endures continuity up to the 2nd 
derivative of the potential. We choose the number fractions of th e two types of particles as x+ = a;_ ~ 1/2 and the 
choose the ratio of the diameters as a+/a^ = 1.2 [3a, IIOTI . Il08l |. The radii cr+ and (t_ are fixed by requiring that 
the effective diameter to be equal to the unit length scale Iq, i. e. (X]e=± e'=± XeXc'{(J^ + a^')'^)^/^ — Iq. Then the 
thermodynamic properties are determined by a single parameter F = plQ^kBT/uo)^^^^ ■ Here p = N/V is the number 
density where V is the volume of the system measured in the unit of Iq . In the following we choose p = 1 and study 
the temperature range 0.05 < T < 0.3 where the temperature is measured in the unit of uo/kB- The mode coupling 
critical temperature Tc and Kauzma nn transi tion temperature Tk of the system has been found to be T^ ^ 0.19 — 0.22 
by molecular dynamics simulations 107 109| and Tk — 0.125 by the cloned liquid theory |19| . 



2. Analysis 

Our analysis proceeds as follows. First we obtain some basic thermodynamic parameters. To this end we analyzed 
the thermodynamics of the system following the scheme by Coluzzi et. al. in [19|. More precisely we used the 1st order 
cage expansion for the binary system which is a straight forward extension of the prescription |18l | summarized in sec. 
IIIC 51 We employed the binary Hyper Netted Chain (HNC) approximation [1^, which is the HNC approximation 
|110| adapted to the binary mixture, to obtain the radial distribution functions g'^'^ (r) between particles of type e = ± 
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FIG. 8. Rigidity of a binary soft-sphere system, computed by the cloned hquid approach. The Ka uzmann t emperature of the 
system is Tk ^ 0.12 [19t | and the mode coupling critical temperature is around Tc ~ 0.19 — 0.22 [l07l ] [l09l |. The intra-state 
rigidity fi (squares) is nearly independent of temperature below Tk while it sharply decreases with increasing temperature 
above Tk. On the other hand, the instantaneous rigidity, i. e. the Born term (filled circles) increases with the temperature. 
The renormalization of the Born term is very weak (open circles). The significant part of the non-afline correction is due to 
the Jl term (triangles), whose magnitude sharply increases with temperature above Tk. 



and e' = ± of the liquid state of the binary system, and the free-energy of the hquid state Eq. (|98p which are needed 
in the thermodynamic analysis. As the main outcome of the thermodynamic analysis, we obtain radial distribution 
functions g^^ (r) , the the glass-transition line m = m* (T) and the cage sizes A^ (T) of the two kinds of particles e = ± 
which are indispensable in the following analysis. 

Next we compute the intra-state rigidity jl of the system at the level of 1st order cage expansion. To this end we 
need to evaluate various terms in Eq. (|152[) by representing them using particle distribution functions as explained 
in Appendix lEl IF] and IGI The values of the parameters c^ Eq. (|G4p . which controls the strength of the non-afhne 
corrections, are obtained using the cage sizes A^^ and the parameter m*{T) determined in the thermodynamic analysis. 
We use the radial distribution functions p*^*^ (r) in the evaluation of the Born term Eq. (|G2p and the non-affine 
correction term J2 Eq. (jG3p . All integrals J dr . . . are done numerically. To evaluate the non-affine correction term 
Jl, which involves three-point particle distribution function (see Eq. (jE12p ). we use the Kirkwood superposition 
approximation Eq. (jE5[) which approximates the three-point function as a product of two-point functions g"^*^ (r) (See 
Eq. (JGIOP and Eq. (|Glip ). To evaluate the non-affine correction term J3 Eq. (|G5|) . we need to evaluate the derivatives 



%r {r;S+,S^)/dS, 



s=o 



In practice, we obtain them as follows: first we evaluate the radial distribution functions 



g'^ ' (r) using the binary HNC [19| for the modified potential Eq. (|G6| with some sufficiently small values of the 
perturbation Se and then we take numerical derivatives of the results with respect to 6^- 



3. Results 



As shown in Fig. [S] the intra-state shear-modulus fi exhibits the following characteristic temperature dependence. 
On one hand it is almost independent of the temperature below the Kauzmann transition temperature Tk ''^ 0.12 
and takes a value around {a, ^ Q. On the other hand it becomes strongly dependent on the temperature at T > Tk. 
Evidently there is a significant difference between the Born term, which represents the affine part of the response, 
and the intra-state rigidity jl at all temperatures. Even at very low temperatures the reduction of the rigidity due to 
non-affine corrections is significant, in agreement with numerical studies of static response to shear of energy minima 
(inherent structures) J62l l63j . Moreover the sharp decrease of (1 above Tk is clearly due to the non-affine correction. 
The fact that the Born term h increases with increasing temperature would appear somewhat puzzling but it simply 
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FIG. 9. Two temperature dependent factors in the non-afBne correction term Ji. In the binary system the Ji term Eq. (|155p 
can be written as (see Eq. (|G7|) ). mJi — X]e=± Ctl3*'S,^. The panel a) shows the inverse of the effective Hookian spring constants 



K^ = l3*Cc with Ce as defined in Eq. (|G4p . multiphed by /3* — m*{T)p. The panel b) shows the squared amplitudes of the H 
fields, HJ as defined in Eq. (|G8|) . 



means that particles, which are interacting with each other by purely repulsive interactions, are colliding each other 
more often at higher temperatures so that instantaneous increase of the energy (without any relaxation) with a given 
shear-strain 7 becomes larger at higher temperatures. 

Let us examine the non-afRne corrections in more detail. As discussed in sec. IIII C 31 the non-affine corrections can 
be divided into two groups: 1) the group of terms ( J2 and J3) which yields the renormalized born term 6°^ Eq. p54p 
together with the original born term 6 and 2) the Ji term Eq. (|155p which is distinct from others. In Fig. |S] we 
included the curves of both the renormalized born term h and the Ji term. Apparently the major contribution 
to the non-afRne correction is the Ji term. As can be seen in Fig. [8l it only weakly depends on the temperature at 
T <Tk but sharply increases with the temperature above Tk- 

From the definition Eq. pSSp (more precisely Eq. (jG7[) which is adapted for the binary case), it can be seen that 
there are two different sources of the temperature dependence of the Ji term. One is the factors k~^ = Ce/3* which 
are the inverse of the effective Hookian spring constants (see Eq. (|147p and Eq. (jlOSp ) and the other is the amplitudes 
of the S fields. As shown in Fig. [HI the main source of the strong temperature dependence of the non-affine correction 
term above Tk is clearly the amplitude of the S field. We have argued in sec. Ill A 81 and sec. IIIIC3l that the non-affine 
correction term can explicitly depend on the temperature because it reflects thermal fluctuations within metastable 
states (metabasins). 



4- Discussions 



We found the intra-state shear-modulus jl becomes almost independent of the temperature below Tk suggesting 
that thermal fluctuations are significantly suppressed at T < Tk- Then the value of fi. '^ 6 at T < Tk may be 
compared with the rigidity of inher ent s tructures ms (see Eq. ([50]) ). Quite interestingly an analysis of ^is Eq. ([57)1 of 
the same system indicates fiis ~ 5 
the onset value of the plateau Goo 
of the same system indicates Goo ^ 



lll| independently of the temperatures. Similarly it may also be compared with 



in the stress relaxation (see Fig. [2]). Indeed a molecular dynamic (MD) simulation 
■^ 5 [70| independently of the temperatures. 

On the other hand, we found the intra-state shear-modulus fi becomes sm aller at higher temperatures at T > Tk 
(See Fig. \TU^. This is very different from both ^is |lll| and Goo [63, [ZO, Ill2| which appear to depend little on 
the temperature up to surprisingly much higher temperatures. Thus we expect Gp ^ /t lies below Goo at T > Tk . 
However, somewhat puzzlingly, in numerical observations of stress-stress auto-correlation function Ca{T) |69l. l70l . Ill2l | 
there are no indications of the existence of such an rigidity Gp lying below Goo- Probably this means that the 
distinction between the /3-relaxation and the a-relaxation is just a smooth crossover in the real dynamics at variance 
with the idealized picture shown in Fig. [2]) . 

In Fig. [ini we compare the intra-state rigidity fi with the static rigidity obtained by a previous MD simulation on 
exactly the same system done by Barrat et. al. (1988) [3al- In the MD simulation [Sg] the static fluctuation formula 
for the elastic modulus (see sec. Ill A2[) was used to compute the elastic moduli. It appears that the theoretical 
curve, which is obtained with no adjustable parameters, compares reasonably well with the data points of the MD 
simulation. However we should note that this comparison is delicate since static rigidity of a truly equilibrated system 
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FIG. 10. Comparison of the rigidity of a binary soft-sphere system computed by Barrat et. al. (1988) [3b | with N — 3456 
and the theoretical one. Here G and G' represent respectively the rhombohedral and tetragonal shear-modulus which become 
identical in isotropic elastic systems. 



of large enough system sizes should become at least down to Tk ^ 0.12 (corresponding to {p(j^ff)^{£/kBT) ^ 8.3 in 
Fig. [TO]). 

It is interesting to note that the intra-state shear-modulus /i, as well as the rigidity observed by the MD simulation 
[Sq . crosses at around T ~ 0.23 (corresponding to (pcrj^ff)^(e/fcBT) ~ 4.3 in Fig. [T0|) , whi ch is rath er close to the 
putative mode coupling critical temperature Tc ^ 0.19 — 0.22 estimated by MD simulations J107f|l09l |. This feature 
may be interpreted in favor of the Born's original criteria for melting |37| , i. e. continuous vanishing of the rigidity 
approaching the spinodal temperature is valid for the present system. The situation is similar to the numerical 
observations of the rigidity by some MD simulations on superheated crystals [391, |40l . However we note that our 
computation is limited to the 1st order cage expansion and does not tell us the location of the putative spinodal 
temperature where the glassy metastable states actually disappear. 

The discussion in sec lIII El suggests rather that the system would actually exhibit discontinuous melting but preceded 
by a continuous, square-root behavior of both the cage size A and the intra-state rigidity /i at a somewhat lower 
temperature. Then the seemingly continuous vanishing of the shear-modulus /i around T^ is a pure coincidence and a 
small discontinuity may be hidden in the numerical data of |36| shown in Fig. 1101 We have discussed two qualitatively 
different driving mechanisms for such a behavior. The strong temperature dependence of the non-affine correction 
term Ji shown in Fig. [5] and the fairly good comparison of jl with the data of the MD simulation shown in Fig. [TU] 
suggest that the coupling between the cage size and the shear-strain fluctuation may be important in the present 
system. In this respect it is interesting to note that the conventional mode coupling theory which does not take into 
acco unt the coupling of the shear-strain fluctuation and the density predicts a too large critical temperature Tc ^ 0.33 
[113|. However we must remind ourselves once more that melting of amorphous solids can be at most just a smooth 
crossover in finite dimensional systems. 



IV. DISCUSSIONS 



Before concluding this paper, let us present discussions on some problems. 

• In the present paper we analyzed static rigidity of the amorphous metastable states of structural glasses. We 
discussed a possible interpretation of the hierarchy of the static rigidities in the context of rheology in sec. Ill A 51 
- sec. Ill A 81 In particular, we argued that the rigidity of the inherent structures pis and the intra-state rigidity 
p correspond respectively to the rigidity at the onset of the plateau regime Goo and the large time limit of the 
plateau regime Gp (see Fig. [5|) . However the plateaus become obscure at higher temperatures and it becomes 
difflcult to separate the /? and a-relaxations. A possible way out of this difficulty is to perform the analysis 
of fluctuation-dissipation relation (FDR) discussed in Appendix IA2I [7C| which allows one to "eliminate times" 
and focus on the breaking point of the FDR, which is expected to correspond to the crossover point between 
the /3 and a relaxations, i. e. Gp. 



• 



• 
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From a broader perspective, the intra-state rigidity fi is just one of many intra-state susceptibilities which 
characterize quasi- equilibrium material properties of amorphous solids, which are well accessible to experimental 
observations and important in practical use of amorphous solids. It is interesting to note that properties of 
the underlying "ground state" of a supercooled liquid can show up in the intra-state susceptibilities even if the 
system as a whole is out of equilibrium. As we discussed in sec. Ill C 4l and in sec. Ill C Ti the cloned liquid approach 
realises decomposition of a generic susceptibility into intra-state (/3-relaxation) and inter-state (a- relaxation). 
This suggests possibilities to make further progress in understanding material properties of amorphous solids 
out of equilibrium starting from microscopic Hamiltonians. 

A natural step to go beyond the mean-field theoretical level is to consider spatial fluctuations of the order 
parameter i. e. the cage size A{v), shear-strain 7(r) and e.t.c. In particular we expect that the coupling of the 
shear-strain 7(r) and the cage size A{v) becomes important in finite dimensional systems as we discussed in sec. 
nil Dl and sec. IIII El Following the discussion in sec. IIII El we may write a phenomenological Ginzburg-Landau 
like free-energy functional of the cloned liquid as, 



• 



F[A(r),7(r)]= d^r 



E '-^^ f ^- - h) ^" w^'(^) + 1 (^^(^))' + ■9™(^(^)) 



(174) 



where p = N /V is the number density, gm{A) = p{Gm{A)/N) is the local 'mean- field' free-energy with Gra{A) 
given in Eq. (|16ip and jJ.{A) is the intra-state rigidity. This can be compared with the free-energy functional 
Eq. (IC9[) of the replica field theory of vectorial spin-glasses discussed in Appendix [C] where the equivalent of 
the shear-strain field 7(r) is the derivative of the angular field V6'(r) which parametrizes the spin-waves. The 
free-energy functional defines an "effective Debye model" for amorphous solids. To be completed the elastic 
free-energy associated with normal strain (compression) must also be added. In principle all key parameters 
which parametrizes the free-energy functional can be determined by microscopic computations using the mean- 
field theory (cloned liquid theory) . Especially computation of the bulk modulus which parametrizes the elastic 
free-energy associated with the normal strains can be done in the same way as done for the shear-modulus fi 
in the present paper. However, we must note that such a coarse-grained description cannot describe important 
now- Debye features in amorphous solids such as the boson peak and anomalous soft-modes around jamming 
point [43. 

We found the rigidity depends on the cage size A as Eq. (|156p which reads as, 

fi{A)^co-ciA + 0{A^). (175) 

This implies the local rigidity /i(r) is smaller in the region wher e the local cage size A{v) is larger. It might be 



related to the numerical observation of the elastic heterogeneity 114l | . It will be very interesting to study more 
in detail the consequences of the coupling between the transverse fluctuation (shear-strain) and longitudinal 
fluctuation (cage size) in the stat ic and dy namic properties of amorphous solids and supercooled liquids, such 
as the dynamic heterogeneity [6^, Ill5l4ll8{ . 

We argued that inter-state part of the response to quasi-static shear brings about jerky, intermittent stress- 
strain curves at mesoscopic scales (see sec. IIIIBI) . They are signatures of the inter-state fluctuations which 
are Goldstone modes in glasses which try to restore the brok en the translational symmetry. It will be very 
interesting to perform some micro-rheological experiment [119l | of slow strain rate in supercooled liquids to look 
for some signatures of the putative quasi-static intermittency. A possible gedanken experiment may be to slowly 
drag a single particle immersed in a deeply supercooled liquid through a Hookian spring attached it. Such an 
experiment c an be realized for instance by applying magnetic field to a magnetic particle immersed in a colloidal 
glass system [ll9[. In order to suppress melting, the system may be confined in a narrow region between two 
walls with separation L comparable or smaller than the mosaic length ^mosaic \A, ISl ■ Then the interior of the 
local region (mosaic) around the particle will behave as a bulk system in the glass phase T <Tk- By dragging 
the end of the spring slowly over some a small distance x, this region will be subjected to a local shear-strain 
of order S"f ex cc/^mosaic- Then the force curve F(x) of the spring will exhibit intermittent profiles either due to 
level crossings amon g di fferent metastable states within the mosaic or yielding due to saddle-node bifurcation at 
metastability limits [33,|33. A useful observable may be the force-force correlation function Cf{x) = F{0)F{x) 
(where "^^^ means averaging over the space or trajectories). Presumably the correlation function Cpix) exhibit a 
cusp like singularity at x = (but rounded over some width Sx^ ex T at finite temperatures) and decays rapidly 
as function of x/Sxg where Sxg is the typical spacing between the shocks. The conjecture Eq. (|133p implies 

OXs ~ ?mosaic07s '^ Smosaic- 
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V. CONCLUSIONS 

In this paper we developed a first principle scheme to evaluate the rigidity of supercooled liquids and glasses by 
analyzing the rigidity of a cloned liquid, which might appear absurd because the rigidity of a liquid is just zero. We 
showed that the replica technique allows one to switch-on or off contributions of the inter-state fluctuations of the shear- 
stresses to find non-zero rigidity of amorphous metastable solids and inter-state responses with intermittent nature. 
We discussed physical interpretation of the result within the context of linear rheology and the phcnomenological 
free-energy landscape picture. We also presented a model computation of the rigidity of the binary soft-spheres and 
found that the result compares well with the result of a previous molecular dynamic simulation. A strong advantage of 
the cloned liquid approach is that it fully makes use of conventional liquid theories such that one can choose the best 
one for a given problem at hand. Since the rigidity is a basic property which is very important both to understand 
the physics of glassy systems and to use amorphous materials in practice, it will be certainly interesting to develop 
the method further and apply it in various glasses and jammed matters. 
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Appendix A: Fluctuation formula for stress relaxation 

1. In equilibrium 

Here we outline a simple derivation of the generic fluctuation formula for the stress relaxation Eq. ([TT|) |54| . To be 
specific we consider the time evolution of the system in terms of Liouville equation of motion, 

||P(t))=£|P(t))) (Al) 

where vector \P{t))) represents probability distribution of the positions r,; and momentum p^ of particles (i = 
1, 2, . . . , A^) and £ is the Liouville operator; 

C = f\ (-^ . A + ^ . AS] . (A2) 

The equilibrium distribution, which is an eigen vector of the Liouville operator with zero eigen value £|Poq)) = 0, is 
given by, 

|Peq))=^-'e-^^"^^e-^^ (A3) 

where Z is the partition function Eq. ([5]). 

The expectation value of the stress a observed at time t can be written formally as, 

(a(t)) = ((<7e^*|P(0))) (A4) 

where (( on the r. h. s of the above equation represents traces over the positions r^ and momentum p.; of the particles. 
The stress is defined in Eq. ([7]) which reads as, 



1 dU 



(A5) 

7=0 



where U is the potential part of the Hamiltonian of the system defined in Eq. ^. We assume that U depends only 
on the positions r^ of the particles. 
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The linear response of the stress with respect to infinitesimal changes of the strain S^{t) can be obtained formally 
as, 

{5a{t)) ^ ((^e"|P(0)M7(t) +£dt'((ae^(*-*')^e^*'|P(0)))<57(t') (A6) 

= {b{t))Sj{t) + J^ dt'{{ae^^'-''^f^e^''\Pm)Sj{t') (A7) 



where 

, da 1 d^U 



d-f N d-y^ 



(A8) 

7=0 



is the Born term defined in Eq. (|137p which represents the instantaneous, affine response to shear. The 2nd term on 
the r.h.s of Eq. (|A7[) represents non-affine responses to shear due to stress relaxations. 

If the system is equilibrated at the beginning, i. e. \P{0))) ~ \Pcq)) we find the fluctuation formula Eq. (fTTj) . 

{5a{t)) = (6)<57(t) -pj^ dt'^^5^{t') (A9) 

where C„{t^t') is the auto-correlation function of the stress, 

C,{t, t') = (a(t)a(t')) = ((ae^(*-*' Ve^*' |P(0))) (AlO) 

Derivation of Eq. (jA9| from Eq. ()A7|) is a standard one: one just need to notice. 



=^^*"*'^^^^*>cc,» =/3((...e^(*-*') [cf. f-CJ e^*>c,» = -/?|,((. . . e^(*-*')^e^*>e,)). (All) 



Here the 2nd equation follows from the fact that |Pcq)) is the equilibrium distribution given by Eq. (jASp . 

2. Out-of-equilibrium 

In out-of equilibrium glasses, i. e. glasses under aging or driven by external forces, the equilibrium fluctuation 
formula Eq. ()A9P cannot be expected. However the linear response can still be written formally as, 

/■* BC (f f') 

{5a{t)} = {b}Sj{t) - 13 j^ dt'X{t,t')^^Sj{t') (A12) 

by introducing an unknown factor X(t,t') which is 1 in the case of equilibrium response. 

As suggested by the studies of the dynamical mean- field theory of spin-glasses [72| , a natural way to investigate 
large time behaviors of out-of-equilibrium glassy systems is to consider large time limits i, i' — >■ c» such that the value 
of the correlation function is fixed to a certain value Ca-(t,t') = C. In such a limit the X factor becomes parametrized 
solely by the value of the correlation function C, 

lim X{t,t') = X{C). (A13) 

Then in the case of the typical stress relaxation protocol discussed in sec. Ill A4[ the linear response in the large time 
limits i, t' — > oo with fixed C = Ca{t, t') can be written simply as, 

(Sad)) i'lin-H^^C^{t,t) 

lim ^—^ = hm {bit)) - 13 / dCXiO (A14) 

The simplest ansatz for the function X{C) would be to assume that it consists of two pieces, 

''^''^^\. {PC <,,/.) ^^^^) 

This is exactly the case for the class of mean-field spin-glass (MFSG) models which exhibit 1 step replica symmetry 
breaking (RSB) including the p-spin spherical MFSG model [22j90]. Because of the intimate analogy between the 
MFSG models which exhibit IRSB and glass phenomenology [a-0, Q: it is often assumed that X{C) fu nction of 
structural glasses also has the form Eq. (jA15|) . This ansatz is supported by numerical simulations J120l4l22 |. 
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Appendix B: Cloned liquid computation of the linear response of p-spin spherical mean- field spin-glass model 

Here let us analyze static magnetic linear response of the p-spin spherical mean-field spin-glass (MFSG) model 
using the cloned liquid approach. The validity of the cloned liquid approach for this system was demonstrated in 
[123| . Here we demonstrate that the decomposition of the intra-state and inter state linear magnetic susceptibility 
can be done exactly as discussed in sec. IH C 41 and also b y usin g the cage expansion approach discussed in sec. IH C 71 



The Hamiltonian of the p-spin spherical MFSG model jl24l | is given by, 



N 
l<ii<i2<---ip<A' i=l 

where the scalar spin variables Si at sites i = 1, 2, . . . , A^ Here we have introduced a uniform probing external magnetic 
field h. are forced to satisfy the spherical constraint X]i=i sf = N. The random couplings Ji^i2...i are statistically 
independent from each other and drawn from a Gaussian distribution with zero mean and variance J'^p\/{2NP~^). In 
the case of p — > oo limit it reduced to the random energy model [92| . 

Let us consider a cloned system with m replicas a = 1,2, . . . ,m with fixed overlap, 



1 ^ 



(B2) 



1=1 



Static properties can be computed from the partition function of the cloned system, which is obtained as. 



a, 6 

where 

32 t2 



GiQab, Xab) = -^ E Qafc - E ^-^Qab + ^ log(dct (2Aa,)/^™) (B4) 

ab ab 

with 

Qab = (1 - q)Sab + q- (B5) 

and 



ab 



In Eq. (jB3P we have simply took average of the partition function over disorder of the random couplings. This is the 
so called as annealed average which is valid only in the high temperature phase where replica symmetry is not broken. 
We analyze the low temperature region choosing small enough m so that we remain in the high temperature phase 
(See Fig. E] a)). 

In Eq. (|B3p wc have introduced different probing field h"' on different replicas. They are assumed to be infinitcsimally 
small. Performing the integration over Xab by the saddle point method { at h"" = 0) wc find, 

lini -/3m ^"'^^;^°^ = -GiQab, Kb) + ^ E ^''Q-"''' ^^b = ^Qab (B7) 

ab 

This is formally the same as the 1 step replica symmetry breaking (RSB) solution |124| . Thus by performing extrem- 
ization with respect to the parameters q and m one recovers the exact IRSB free-energy of the system [l23J |. 
Using the above results the linear-susceptibility matrix Xab Eq. (j80p is readily obtained as 



Xab = PQab = /3(1 - q)6ab + Pq. (B8) 

Comparing with Eq. (j80p we can easily identify the intra-state linear susceptibility x = /3(1 — q) and inter-state 
susceptibility^ = mjiq. The total linear magnetic susceptibility becomes 

m 
b=l 
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As expected these results agree with those of the usual 1 step RSB solution [l24l | . In the standard interpretation of 
the RSB solution l9l| for the spin-glass phase below the static spin-glass transition temperature TgG: the factor m 
(written rather as x in the context of RSB in quenched disordered systems) in front of the 2nd term on the r. h. s 
is understood as the probability that two replicas stay at different metastable states. More physically the factor m 
roughly corresponds to the probability that the system stay at an low lying 'excited' metastable state rather than 
'ground' metastable state in the glass phase. The parameter m is found to behave typically behaves as to ~ T/Tsg 

In the temperature range Tqq < T < T^ (corresponding to Tk < T < Tc), the system in the liquid state in the 
thermodynamic sense. However the RSB solution Eq. (|B9|) is still very useful if we consider m —^ 1 limit in the following 
sense. First note that the total susceptibility becomes lim„j_j,]^- x = /? so that the correct paramagnetic susceptibility, 
which could be obtained by the replica symmetric (RS) ansatz which yields q = 0, is recovered independently of the 
value of q. The value q should be interpreted as self-overlap of metastable states which can be obtained as lim,„_j.]^- q 
using the RSB ansatz or equivalently the cloning liquid like methods [23, [Sj, Il23{ . Indeed the latter procedure 
correctly reproduce the self-overlap value q of metastable states which dominates the thermodynamics of the "liquid 
state" at Tsg < T < Td- Thus we can still decompose the total susceptibility into the intra-state susceptibility 
X = /3(1 — \i'aijy^__^i- q) and the inter-state susceptibility x = /^limm-s-i- Q in this intermediate temperature range. 

Now let us try the cage expansion approach to evaluate the linear susceptibilities as discussed in sec. Ill C 71 To 
this end we decompose the spin variables into two parts as. 



where Si is the 'center of mass' of the spin 



s^ = S., + 5s1 (BIO) 



m 

TJ ^ ^ 



m 

a=l 



and (5s" is the molecular coordinate which describes fluctuations around the 'center of mass'. As far as the spherical 
model is concerned, the effective action is already quadratic in terms of the spin variables (see Eq. (jB6[) ) so that the 
decomposition Eq. (jBlip is not particularly needed in order to compute the free-energy. Nonetheless, examination of 
the linear susceptibility by this approach may be instructive as we outline below. 

The Hamiltonian of the cloned system can be written in terms of the new variables as, 

m 
l<ii<i2<...ip<N a—1 i a i 

where we have introduced probing fields conjugated to the 'center of mass coordinate' 5'^ and 'molecular coordinate' 
Ssf. Then we only need to replace Eq. (|B6p ) with, 

V-"Vdct(2Aab) J \^ ^^ 
= 7r + t;} ^Sh [l-mdab) oh (B13) 



2 2 "^-^ m 

ah 



Now we can readily find, 



(Js'^Js") = -(1 - m5ab)^-^. (B14) 



This result is reminiscent of Eq. ([92|) combined with Eq. (|102p ) . We can notice indeed that 1 — g plays precisely the 
role of cage size A. Using the above results the linear-susceptibility matrix Xab can be obtained as, 

Xab = iSis^s") = /3{{s + Ss''){s + 5s'')) = f3{S^) + (3{5s''Ss'') = /?[(1 - q)5ab + q]. (B15) 

The last result agrees with Eq. (jBSP obtained without using the cage expansion. 



Appendix C: Rigidity of spin- waves in vectorial spin-glass systems 

Here we present a simple mean-field theoretical analysis of the rigidity of spin- waves in a class of vectorial spin-glass 
systems. Our purpose is to extend the analysis of the rigidity of spin- waves in the superheated ferromagnets discussed 
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in sec. IIIB II to the case of some spin-glass systems approaching the dynamical tran sition temperature Td- To this 
end we consider a replica field theory with rotational invariance following Ref. 104j which analyzed the case of the 
Sherrington-Kirkpatrick model with vectorial spins. 

Although the following discussion does not depend on the details of specific choices of microscopic models but we 
may consider as a specific example the p-spin spherical MFSG model Eq. (JB1[) extended to a model with vectorial 
spins, 



l<ii<i2<...ij,<N 



•-'iii2...ip[Sii ' Si2)(Si^ 



1 j ■ • ■ (Sip-i ■ Sip j. 



(CI) 



The spin variables are two-component vectors Si = {si^i,Si^2) and the symbol • represents the inner product, i. e. 



Z^u=l 



2 •'^^''J.A'- 



Note that the model has the 0{2) symmetry just as the ferromagnetic system considered in 



sec. IIIB II The vectorial mean-field model can be solved by introducing a tcnsorial order parameter, 



1 ^ 

Qab = jj / A^i^si 



b ) 

■i.-ul ■) 



(C2) 



4=1 



where (...) stands for thermal averages. 

This class of systems with p being an even integer greater than 2, i. .e p = 4, 6, 8, . . . exhibit the dynamical transition 
at T d and static spin-glass transition at Tsq at a lower temperature just as the original scalar p-spin MFSG model 
Il24| discussed in appendix [BJ The case of p = 2 is the vectorial version of the Sherrington-Kirkpatrick model studied 
in |l04| which exhibits the usual spin-glass transition at a single critical temperature. 

To analyze generic properties of spin-waves in system as the one given above, let us consider a replica field theory 
of m replicas (a = 1,2,..., m) with a schematic free-energy functional. 



f[q: 



fJ.U-i _ 

b\ ^ 



d'^ 



iEE(vra^ 

a<b liu 




(C3) 



Here the order parameter Q^^(r) is considered as a continuous field in a d-dimensional space. Just as considered in 
appendix [B] we assume that this is a cloned system: replicas are forced to belong to a common equilibrium state. The 

local potential w wZni> {Q^b) ) ^'^ *^® integrand of Eq. (jC3[) can be computed as the mean-field free-energy of the 
underlying microscopic model such as the one given by Eq. (|Cip but the details are not important in the following 
discussion. 

As noted in |104| . an important symmetry property of the free-energy functional Eq. (jC3p is that it is invariant 
under global rotations of any one (or more) of the replicas: it is invariant under a generic transformation of the form. 



Ql 



fj.' u' 



{R 



t\fj.fj.' I 



')a)Q''J i?" 



(C4) 



where 9a and Oh are the angles of the global rotations of the replica a and h respectively. The rotations are represented 
by the rotational matrix i?* {9) , 



R\9) 



_ ( cos{9) ~ ain{9) 
sin(6l) cos(6l) 



(C5) 



Then by choosing a specific reference coordinate system for the 2-component spin space in which the matrix Q^^ 
becomes diagonal, i. .e. Q'^'^ = Qab^ab, a generic replica field Q^^ can be parametrized as. 



Qab = qabJ2^RT''iOa)R'"'{eb) = QabiRT^iOa - 9b). 



(C6) 



Using the latter parametrization the free-energy functional Eq. (|C3[) becomes 



F[Q'ab] = 



d'^J 



a<6 a<6 



(C7) 



Note that the latter two terms of the integrand appear also in the usual scalar replica field theories. On the other 
hand, the 1st term can be naturally interpreted as the spin-wave term just as in the case of the ferromagnetic model 
(See Eq. dMl)). 
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As discussed in appendix [B] the matrix qab of the cloned system can be parametrized as, 

qab = {l-q)6ab + q (C8) 

with q being the self-overlap, i. e. the Edwards- Anderson order parameter. Then the free-energy functional can be 
simplified as, 



F[QZ]^ jd'^r 



a, 6 a<fe 



(C9) 



Here ^ab is the rigidity matrix of the spin-waves which is obtained as. 



^lab = cmq^ [Sab I • (CIO) 

V "V 

Remarkably the rigidity matrix Eq. (jC10[) has precisely the same structure as that for the shear-modulus of glasses 
(see sec. IIII Ap : /i = fiSab + A with jl = —{\/m)jl so that the sum rule YlT=i ^^a■b = is satisfied. Now the 'intra-state 
spin-wave rigidity' can be easily read off as, 

(i = cmq^. (Cll) 

Evidently this is quite similar to the case of ferromagnet discussed in sec. IIIB ll in the sense that the spin- wave rigidity 
is proportional to the square of the order parameter (here the Edwards- Anderson order parameter q) . 

We readily know that the Edwards-Anderson order parameter q jumps discontinuously to approaching the dy- 
namical transition temperature Tj from below, preceded by a square-root singularity. 



q{T) - q{Td) cc ^Td - T. (C12) 

Using the latter in Eq. (jClip . we immediately find that the 'intra-state spin- wave rigidity' should discontinuously 
vanish approaching Tj from below preceded by a square-root singularity. 



/i(r) - A(Td) ex cmy/Td - T, (C13) 

where 771 = 1 at Tsg < T < Td- Apparently these features are the same as the superheated ferromagnet discussed in 
sec. IIIB II 

Appendix D: Proof of the relation Eq. (|148p 

Suppose that the original 2-body interaction potential v{r) (See Eq. (|4])) is slightly deformed as v{r) — >■ v{r) + 5'w{r) 
where u;(r) is an additional potential and (5 is a small parameter. Notice a trivial fact that the rigidity Eq. (|S]) must 
be zero irrespective of the deformation of the potential as long as the system remains in the liquid phase. Then the 
identity, 

{b)~(3[{a^)-{ar~]=0, (Dl) 

must continue to hold under variation of the parameter S. Especially the form of the renormalized potential Vcs{r) 
given in Eq. (p4)) suggests us to consider the case w{r) = V^w(r). 

Let us examine changes of the thermal averages on the 1. h. s. of Eq. (JD1[) under variation of the parameter S. From 
the microscopic definitions of the shear-stress a and the Born term b given in Eq. ([7]) and Eq. (|9]) , we find that they 
are actually derivatives of the potential v{r) so that they change under variation of S. Thus the thermal averages can 
expressed formally as, 

r ^ rid r ^ rid 

lr\ \ '7-1 / TT ''* -;3y"_ u(r)+(5tu(r)^ ry / 77*^ ^* -S Y" .^ . ti(r)-|-5tu(r) ,'-^\0^ 

■'^ i=\ ■'^ i=\ 

with O being the Born term 6, the shear-stress cr or cr^ . Expanding Eq. (|D2|) in power series of b and inserting the 
results into Eq. (jDip . we can find an identity at each order of b. 
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Choosing the additional potential in particular as w(r) — V^v(r), we find at 0{S) the following identity, 

i<j i<j k<l 

-2^ E E /?[(V^^(r,,,Ja(r,,,J)-(VV(r,,,J)(a(r,,,g)] 

il<il i2<J2 

+ ^ E E E/3'^[Hr^u■J^(r..,JV2^;(r,0)-Mrn,,Mr..,J)(V2«(rfe^))] (D3) 

il<il i2<J2 k<l 

Using the above result for the liquid at temperature T/m*, we find the relation Eq. (|148p . 

J2 + J3 + m*{'h + h) = 0. (D4) 

Appendix E: Representation in terms of particle distribution functions 

In order to compute each terms in Eq. (jl52p . we need to evaluate thermal averages (. . .), over the thermal fluctu- 
ations of several quantities which are multi-point functions of the CM coordinates of the molecules. To this end it is 
convenient to represent each terms in Eq. (I152p and related quantities in terms of particle distribution functions, 

(El) 

where p = N jV is the number density. Here we introduced a short-hand notations X^istj ~ Ei Sill ^'^y)' ^Tlii=j^i=tk ~ 
Si X^i X]fc(l ~ SijSjk). Since we will consider homogeneous and isotropic systems, we can write, 

p'9ir)-^Y.{^'ir^j^r)) (E2) 

CXD / 

p'g^'Hr,r') - ^ E (-^'ir., - r)S'{r,, - r')) = J29l'\r,r')P,{cos{e)) cos(0) = ^ (E3) 

i^j^k 1=0 

In the first equation g{r) — g^'^\r) is the usual radial distribution function. In the last equation Pi{x) {I = 0, 1, 2, . . .) 
are the Legendre polynomials {Po{x) = 1, Pi{x) ~ x,P2{x) = ^{3x^ — l)v) and 6 is the angle between vector r and 

r'. The coefficients g^ {r,r') are defined as, 

9^i'\ry) ^ -^ J\9sm{e)Piicos{e))g(^Hr,r'). (E4) 

The radial distribution function g{r) can be comp uted using the standard closures of the liquid theory [23| such as 
the hyper-nettcd-chain approximation (HNC) [llOJ and e.t.c. A simple way to evaluate the three-point correlation 
fmiction ^"'(r, r') is the Kirkwood superposition approximation, 

g(3)(r,r') = ff(r)g(/)ff(|r-r'|). (E5) 

1. Cage size 

The cage size A at the 1st order cage expansion Eq. (|104p of d = 3 dimensional system can be expressed as, 

3 

Pp L dr4:Trr^g^(r)'S/'^v{r) 

with the radial distribution function 5*(r) at temperature T* = T/m*(T) and, 

7,(1) 

V^wfr) :=:w(2)(^)^2^— . (E7) 

r 
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2. AfRne response 

Similarly the Born term can also be expressed just by the radial distribution function g{r) as, 

/^ 1 2 (2) 4 



oo 



drAnr g^{r) 



Vl5 



,(1) 



15 



(E8) 



where we used 6(r) = r^v^'^'x^z^ + rv^^>(l — x^)z^ given in Eq. ([9|). Here (. . .)angic denotes an "angular 
average" over isotropic fluctuations of the orientations of the unit vectors {x, y, z) = {x/r, y/r, z/r) with r = 
y/x^+y^+z^. 

The terms J2 and J3 are corrections to the Born term due to the rcnormalization of the potential Eq. ([M|) . The 
term J2 just involves 2-particle distribution function and we find, 



P 



J, = -cl Y. i^'bir.,))^ = -cf J^ drAnr'g^r) 



i<j 



—r^iM\r) + -rf(3) r + —v^^\r) + — ^ 

15 3 15 15 r 



(E9) 



On the other hand J3 is more complicated since it involves connected 4-particle distribution function. This term 
can be evaluated by considering a liquid with a modified potential. 



Vcii{r] 5) = v{r) + (5V^w(r), 
at temperature T* = T /m. Let us denote its radial distribution function as (7,(r; 5). Then 

^3 = 4EEKKr,,)/3*V\(r,0>* - (&(r.y)>*(rv2t;(rfeO>*] 



(ElO) 



.P l~ j„._ 2 dg*{r]5) 



Jo 



i<j k<l 

o 

drAirr' 



dS 



S=Q ^ Jo 



dg*ir;S) 



dS 



6=0 



15 15 



(Ell) 



3. Non-affine response 

Next let us examine the mJi term defined in Eq. (|146p . By noting J2i Yl,ji{^i) Sj2(#i) ^ ^i^j '^^t^nT^]2' *^® 
term mJi can be decomposed into two parts as, 

"^^i=4EE E /^*(^^('-'^i)-^'^(r^^-^))*=4E/?*(l^^'^-|')*+4 E r(Va(r,,J.Va(r,,J), 






s¥j 



s¥ii#J2 



C/3V / dV5*W|S(r)|^+c/3V2 / d'^nd'^rs^f '(ri, r2)S(ri) • 5(r2) 



(E12) 



rnJii mJi2 

Using Eq. ([7|) wc find the explicit expression of the vector S — (S^;, Sy, S^), 

S.(r) 



^^ - [.,,(2) 

^^ = M2)(r)z2+t,(i)(l-z2)]^. 



[rw^ ^(r) ~v '{r)]xyz 



(E13) 



The 1st term in the r.h.s of the last equation of Eq. (|E12p . i. .e mJn can be easily evaluated using the radial 
distribution function g{r), 



mJii = c/3> / d'^rg,{r)\E{r)\^ = c/3*p / drAirr^g^r) 



-rz;(2)(r))2 + IrtS^) {r)v^^) (r)) + -v^'\r)f 
15 15 15 



(E14) 
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The 2nd term in the r.h.s of the last equation of Eq. ()E12[) . i. .e toJi2 involves the 3-particle distribution function 
gi-^)(^r,r'). As shown in Appendix |F] it can be cast into the following form, 



™Ji2 = c/3*p2y d<i^^d<i^2r?rlgi''(ri,r2)S(ri) •S(r2) 



„2„2 



c/3*p I dridr2r2rir2 
175 



9i )*{ri,r2)—^ ( v^ >{ri)j ( i;^ •'(rs) 

(E15) 



+ iAin,r2)^inv'''\n) - v^'\n))ir2v'''\r2) - v^'\r2)) 



where g\ (ri , r2) and g^ {ri , r2) are the coefficients of of the expansion of the three particle correlation function 
5^^-'(ri,r2) by the Legendre polynomials (see Eq. (|E3[) ). 



Appendix F: Evaluation of mJi2 
Here we show derivation of the 2nd equation of Eq. (jElSp . The term mJi2 is defined as, 

m Ji2 = cl3*p^ J d''nd''r2gi''^ (n, r2)H(ri) • H(r2). (Fl) 

By using the formal expansion of the three particle distribution function by the Legendre polynomials, 

5(3)(ri,r2) =^gp)(ri,r2)P,(cos0) cos0= ^ 



1=0 



and a formal expansion of the field H(r) by spherical harmonics Y^ 



(F2) 



(F3) 



where 9 and are the angular variables of the 3-dimensional polar coordinate r = r(sin0cos0, sinflsini/), cos0), 
Eq. (jFip can be rewritten as. 



^Ji2^cP* j dr, j dr2rlrlY,{gf\{r,,r2),^ ^ CT {v) ■ [C^ {r2)r 



(F4) 



where [...]* means the complex conjugate. 

Explicit expressions of the components of the vector field H() is given in Eq. (jE13| . From the latter we find, the 
coefficients C™ of their expansions by the spherical harmonics as, 

(Cr(r)). = (Cr(r)), =5,i5,„,oyf {i (-™(2) +.(!)) -.(i)| +<5,3'5„,o/f{-^ (-r.^^' +.«) 



^A rA A ^ /4^/%/5! 

V 7 I oO 



-r«(2)+t,(i) 



(Cr(r)), = V^5i,,{5ma + 5™.-2) J y I ^ f-™(2) 



(_™(2)+^(l)U. 



(F5) 



Using the above result in Eq. (|F4p we obtain Eq. (JEISP which reads, 

{9f)*{'ri,r2) 



mJi2 = cj3* p^ i dridr2rlrl 



167r 



(3, ,^ „,32.V .W(n)-r,.(^)(.0 ^,,,^\ f!!^!M^lH!!!!lM _,(!)(,,) 



9 



+ {9i\iri,r2)-^iriv^^Hn) - v^^\n))ir2V^^^ {r2) - v'^'^\r2)) 



175 



(F6) 
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Appendix G: Formulations for the binary mixture 

Following Coluzzi et. al. in [l9j we denote the radial distribution function between particles of types e = +, — and 
e' = +, — as g*^*^ (r). At the order of 1st order cage expansion, the cage sizes of the two types of the particles are 
found to be, 



A, 



Ppj^ drAnr^ Y,,, x,,gf {r)V^v^^' (r) ' 



(Gl) 



which is a generalization of Eq. (|E6p . 

The Born term Eq. ()E8P becomes in the binary case. 



/•OO / -1 A 

Similarly the J2 term Eq. (jE9p becomes, 

/•OO 

J2 = -Vc,^/ dr4TTr^XeXe'gl'''{r) 



15 



where the parameters c^ are related to the cages sizes A^ as. 



lr\v-')(^\r) + ^r(."')(3) (.) + |(."')(^)(r) + ^^^^^^ 



Ce =2- 



A, 



which is a generalization of Eq. (|147p . Similarly the J3 term Eq. (|E11[) becomes. 



J3 = V] Q- / dr47rr^ V] : 
^ "'0 ,,,„ 



agr"(^;^+,^-) 



d5e 



5=0 



1^2( .'e")(2) £ ( .V')(l) 

15 ^ '^ 15 ^ ' 



where g*^ ' (r; (5+, J_) is the radial distribution function of a system with modified potential, 

«3 = + ,- 



The non-afiine correction term Ji Eq. (|155p becomes, 

m Ji = ^ Q/3*S2 



(G2) 



(G3) 



(G4) 



(G5) 



(G6) 



(G7) 



e=± 



where 



-Ni-E{iEv.K 



(G8) 



i6£ \ i(#» 

with iVe = A^Xf. Using the particle distribution functions we find Jn Eq. (|E14[) and J12 Eq. (|E15|) become 
mJii = ^ Q/3*p / drAnr^ ^ x^x^'gl" (r) 

TO-^12 = ^ Ce/3*p^ / dridr2rlrl ^ [ 

e -^ e'e" 

32^2 /(^«')(1)(^^) _ ^^(^«')(2)(^^) 



^^(„«')(2)(,))2 ^ ^,(„ee')(2)(^)„(l)(,)) _,_ ^ (^«' )(1) (,)) = 



(G9) 



((g"'^")f)).(ri,r2): 



(„-)(.(,)) U^-'T\r.)-r.^v-'r\r.) _ (,e.')(r)(,^; 



IGtt^ 



+ ((5"'^'')^'^)*(n,r-2)^^(ri(«"')('n^i)-(«"V'H^i))(r-2(«"'')(')(r2)-(^;"'V'H^2)) 



(GIO) 
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where (g*^*^ '^ )l (ri , r2) with I = 1,2,... are the coefBcients of the expansion of the three-particle distribution function 
by the Legendre polynomials (See Eq. (jESp V Using the Kirkwood superposition approximation Eq. (|E5P we evaluate 
them as, 

(g"''")f)(r,/) = ^y"'desin(0)P,(cos(0))5"'(r)/^"(r')/''(|r-r'|) cos(0) ^ ^. (Gil) 
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